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ATTACHMENT  II 


ABSTRACT 


The  parameters  in  a  linear  regression  model  can  be  estimated 
by  minimizing  the  sum  of  the  absolute  residuals  (L^  estimation)  instead 
of  the  more  classical  approach  of  minimizing  the  sum  of  squared 
residuals  (  least  squares  estimation).  In  addition  to  other  nice 
properties  estimators  are  less  sensitive  to  outliers  than  least 
squares  estimators.  This  paper  describes  a  linear  programming 
algorithm  and  computer  program  for  obtaining  unbiased  estimators 
and  estimates  of  their  covariances.  These  estimated  covariances 
are  the  new  feature  in  this  work  and  are  an  extremely  important 
ingredient  in  hypothesis  tests  and  confidence  interval  construction. 
Technical  Report  65  provides  an  analogous  treatment  of  estimation 
subject  to  linear  constraints  on  the  parameters. 
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Unbiased  Estimators  and  their  Covariances 


1.  An  Introduction  to  MRS.  A. 

Consider  the  linear  regression  model  in  the  form 

y  -  X0  +  e  (1) 

where  y  is  a  vector  of  n  observations,  X  is  an  n  x  p  matrix  of  rank  p  of  ' 
known  constants,  0  is  a  vector  of  p  unknown  parameters  and  e  is  a  vector 
of  independent  random  variables  (noise)  symmetrically  distributed  with 
mean  zero  and  variance  o2  .  Unbiased  estimation  of  0  can  be  obtained  under 
several  different  optimality  criteria.  The  classical  least  squares  approach 


is  to  estimate  0  by 


0  -  (XTX)_1  XTY 


which  has  the  smallest  variance  among  the  class  of  unbiased  linear  functions 

«• 

of  v.  The  least  squares  estimator,  0,  is  extremely  sensitive  to  large 
values  of  |e|,  outliers,  particularly  when  the  sample  size,  n,  is  small 
relative  to  p,  say  n  _<  2(p+l)  .  This  sensitivity  suggests  that  an  optimality 
criteria  other  than  minimum  variance  should  be  considered.  Several  authors 
(Barrodale  (1968),  Charnes  and  Cooper  (1964),  Gentle,  Kennedy  and  Sposito 
(1977),  Harris  (1950),  Harter  (1974),  Rice  and  White  (1964),  Taylor  (1973)) 
have  suggested  that 

n- 

*  I  y.  -  X.0  |  (3) 

1-1  1  1 

should  be  minimized  with  respect  to  0  where  y^  is  the  i-th  observation  and 
la  the  i-th  row  of  X.  The  estimator,  0,  which  minimizes  the  sum  of  the 
absolute  residuals  is  often  called  the  estimator. 

Since  the  estimate  is  not  necessarily  unique,  the  unbiasedness  of 
an  estimator  depends  upon  its  method  of  computation.  Hartley  and  Sielken 
[1973]  have  shown  how  to  obtain  an  unbiased  estimator  using  any  conventional 


i 
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linear  programming  algorithm  and  an  initial  unbiased  antisymmetrical 
estimator  of  8,  say  8q  ,  where 

8  -  80(e)  -  -  (8  -  60(-e)]  .  (4) 

The  computer  program  MRS.  A.  implements  an  algorithm  for  Minimizing 
the  Sum  of  the  Absolute  Residuals.  The  algorithm  uses  the  least  squares 
estimator  B  as  the  initial  unbiased  antisymmetrical  estimator. 

It  is  nice  to  be  able  to  compute  an  estimate.  The  fact  that  the 
estimator  can  be  made  unbiased  is  the  first  step  in  understanding  its 
properties.  The  second  step  is  to  estimate  its  covariance.  Such  an 
estimate  would  usually  be  a  prerequisite  for  confidence  intervals  or 
hypothesis  tests.  In  the  past  the  absence  of  such  a  covariance  estimator 
has  made  estimation  less  attractive.  MRS.  A  contains  a  mini-Monte 
Carlo  procedure  for  estimating  the  covariance  of  the  estimator.  This 
feature  sets  MRS.  A  apart  from  other  estimation  procedures. 
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2.  Computational  Procedure 

The  problem  of  minimizing  the  sum  of  the  absolute  residuals  can  be 
formulated  as  follows: 


□ 

min  E  r  (5) 

i=l 


subject  to 


i  i  s  1|  •••*  n  , 


riiO  , 


0  unrestricted  , 


(6) 

(7) 

(8) 


where  r^  is  the  i-th  absolute  residual.  However,  to  insure  that  the  re¬ 
sulting  estimator  is  unbiased,  the  problem  is  reformulated  following 
Hartley  and  Sielken  [1973],  In  particular.  Introducing  the  antisymmetrical 
least  squares  estimator,  transforms  (6)  to 

~ri  -  yi  -  v„ "  xi 05  ■  v 1  ri  • 1 " 1 . .  (9) 


Then,  using 


8  -  0O>  -  e(2)  , 

-8  -  0(1)  -  0(2) 

0  0  0 


with  0(1),  e(2),  8(1)  ,  e(2)  >  0  in  (9)  yields 

0  0 


-X.  (0O)  +  0<2))  +  X  (0(2)  +  0(1))  -  r  <  -  y  +  X  0 
1  0  X  0  X  X  1  A 


XA  (0(1)  +  0(2))  -  xt  (0<2)  +  0(1>)  -ri  <y±-  XA0 


t 
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or  equivalently 


-X.  B,  +  X  B,  -  r.  <  -  y.  + 

1  ?  1  1  q 


Xi  B1  "  Xi®2  "  ri  -  yi  "  V 


(10) 


for  i  «  1,  ....  n  where 


B-  -  +  B^2^  >  0 

x  ft  — 


B.  «  6(Z)  +  6(1)  >  0 
2  0  ~ 


(ID 


Now,  in  order  to  compensate  for  any  idiosyncrosies  in  the  particular 

I 

linear  programming  algorithm  used  to  solve  the  problem  in  (5),  (7), 
(10),  (11),  the  problem  is  considered  in  two  equivalent  symmetrical 
forms,  P^  and  ,  with  MRS.  A  randomly  selecting  either  P^  or  P^  with 
probability  %  .  The  problems  P^  and  P^  in  matrix  notation  are  as 
follows: 


P. :  min  E  r . 

i-1  1 

subject  to 


-  - 

—  — 

r  ~ 

-X  X  -I 

B1 

-y  +  XBo 

< 

X  -X  -I 

B2 

_ y '  xv 

B1  »  B2  »  r  —  0  * 


P«:  min  E  r 

i-1  3 


subject  to 
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- 1 

1 

X 

X 

1 

B2 

< 

L  x  -x  -i 

B1 

-y  +  XBo  _ 

r 

Bi  ,  B2  ,  r  0. 

After  either  P^  or  P2  has  been  solved  by  MRS.  A  ,  the  unbiased 
estimator  is 

3  *=  B.  -  B,  +  0  .  (12) 

1  2  o 

2 

If  the  sample  size  was  quite  large,  say  n  >>  (p  +  1)  ,  the  sample 

A 

could  be  randomly  subdivided  into  G  groups,  then  B  estimated  for  each 

A 

group  g  separately,  and  the  covariance  of  B  estimated  from  the  sample 

covariance  of  the  3  's. 

g 

.  MRS.  A  estimates  the  covariance  of  3  using  a  mini-Monte  Carlo 
procedure.  Conceptually  a  Monte-Carlo  estimate  would  be  obtained  by 
generating  several  sets  of  n  y's  ,  finding  the  estimate  for  each 
set,  and  computing  a  sample  covariance.  There  are  two  difficulties 
with  generating  the  y's;  namely, 

(1)  3  is'  unknown,  and 

(2)  c2  is  unknown. 
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The  first  of  these  difficulties  can  be  overcome  by  expressing  the  sum  of 

absolute  residuals  as 

n 
£ 

i-1 

n 

-  £  I  X  B  +  £  -  X  e  I 

i-1  1 


-  £  I  e.  -  X  (6  -  6)  I 
i=l  1 


-  £  |  e  -  X  | 

i-1 


(13) 


using  (1)  and  60  =  0  -  0  and  noting  that  the  covariance  of  the 

A 

estimator  of  60  is  the  same  covariance  of  0  .  Thus,  only  sets  of  n 
e's  need  be  generated.  To  deal  with  the  second  difficulty  (the  unknown 
o2),  note  that 

.  “ 

£  |  e  -  X  60  | 

i-1  X 

-  o  £  |  e*  -  X.  60*  |  (14) 

i-1  1 


*  A 

where  60  -  (0  -  0)  /o  and  /  a  is  symmetrically  distributed  with 

mean  0  and  variance  1.  Furthermore, 


Covariance  (0)  -  Covariance  (60) 

2  ,  *  *  , 

-  ^Covariance  (60  ). 


(15) 


Thus,  the  mini-Monte  Carlo  procedure  generates  K  sets  of  n  e  's;  finds 

A  ^  A  A 

the  K  estimates  6$  ;  and  then  estimates  the  covariance  of  B  by  o 


(9 
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tlmes  the  sample  covariance  of  the  6(3 


Several  possibilities  for  o2  have  been  considered 


(1)  o 


2  n 


-  E  (y  -  X  5 
1  i=l  1  1  C 


)  /  (n-p) , 


(2)  o  2  =  E  (y.  -xiw  (n-p), 
2  i-1  1  1 


(3)  a  2  -  E  (y  -  X  0  )2  /  (n-p) 
3  1  1 


r*  (,*k  -x±  sK>2 '  <-»• 

k=l 


(4)  o  -  E  |  y±  -  X16  |  /  (n-p) 

-  _ _ . 

k .  A  i  4  -  *i  sK  i  '  <»-">  | 

k-l  1=1 


Of  course,  o  Is  the  usual  least  squares  estimator  of  a  while  a 
1  2 

has  the  same  form  but  uses  the  estimator  of  (3  instead  of  the  least 

-  2  ■»  2 

squares  estimator.  The  ratio  estimators,  o  and  o  reflect  the  fact  that 

3  4 

2  * 

the  variance  of  the  yrs  is  o  while  the  variance  of  the  e  *s  is  1  • 

Since  estimation  is  used  to  avoid  certain  weaknesses  in  least  squares 

2 

estimation,  it  seems  appropriate  to  estimate  a  using  the  absolute 
*  2 

residuals,  as  in  o  ,  instead  of  the  squared  residuals.  Furthermore 
4 

E  [|y  -  X6 | ]  is  often  proportional  to  o.  In  particular. 


E  { |y  -  Xg|]  -  o  (2//T)"  , 


If  e  has  a  uniform  distribution; 
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-1 

E  C |  y -  XB| I  -  a  (.472)  , 


if  c  has  a  normal  distribution;  and 


-1 

e  1 1  y-  XB|]  =  a  (/T)  , 


(17) 


(18) 


if  e  has  a  double  exponential  distribution.  Of  course,  if 

E  [|  y-  XB[1  =  Co  (19) 

and  the  e  's  are  generated  from  the  same  distribution  but  with  variance  1, 

-  2 

then  the  proportionality  constant  doesn't  affect  a  .  For  these 

4 

*  2 

reasons  MRS.  A  uses  o 

4 


The  results  in  (16),  (17),  and  (18)  suggest  another  alternative 

A 

for  a  ;  namely, 

o  -  7  E  |  y  -  X  8  |  /  (n-P) 

5  °  i«i  1  1 

•  A 

where  C  depends  on  the  assumed  distribution  of  e  .  Note  that  o 

5 

A 

reflects  only  the  variability  in  Y,  -  X.8  whereas  o  reflects  not 

i  i  4 

only  this  variability  but  also  the  variability  associated  with  the 
linear  programming  algorithm.  Some  empirical  behavior  of  these  five 
estimators  of  o  is  reported  in  Table  1. 

MRS.  A  allows  the  user  to  generate  the  e  s  from  either  the 
uniform,  normal,  or  double  exponential  distributions.  These  distributions 
were  selected  as  being  representative  of  short,  medium,  and  long  tailed 
distributions  respectively.  These  three  distributions  are  also  interesting 
because  maximum  likelihood  corresponds  to  minimizing  the  maximum  absolute 


Table  1 


Empirical  Results  on  Five  Estimators  of  o 


For  n  1 

o2 

=  20 

**• 

H 

*2 

°4 

-2 

°5 

e  distribution 

25 

21.14 

21.66 

20.36 

24.24 

23.45 

Uniform 

25 

21.59 

21.66 

20.26 

22.58 

23.61 

Normal 

25 

40.53 

45.47 

48.77 

36.22 

36.50 

Double  Exponential 

100 

111.63 

119.44 

109.54 

110.12 

113.59 

Uniform 

100 

87.64 

87.92 

81.21 

91.21 

95.67 

Normal 

100 

38.25 

40.06 

34.43 

39.20 

45.93 

Double  Exponential 

400 

375.31 

376.33 

351.39 

360.95 

362.25 

Uniform 

400 

349.55 

350.70 

324.38 

363.56 

381.67 

Normal 

400 

333.70 

336.78 

300.91 

214.12 

251.63 

Double  Exponential 

2500 

3255.33 

3782.59 

3551.58 

3762.09 

3724.80 

Uniform 

2500 

• 

2174.03 

2181.34 

2029.00 

2270.24 

2378.90 

Normal 

2500 

2480.54 

2907.44 

2497.52 

2345.71 

2720.68 

Double  Exponential 

8100 

8363.82 

9644.51 

9027.74 

8310.63 

8470.63 

Uniform 

8100 

7117.80 

7139.82 

6681.97 

7433.47 

7765.47 

Normal 

8100 

7987.26 

8063.69 

7608.09 

8355.12 

9082.50 

Double  Exponential 
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residual  if  the  e’s  are  uniform,  minimizing  the  sum  of  squared  residuals 
if  the  e’s  are  normal,  and  minimizing  the  sura  of  absolute  residuals  if 
e’s  are  double  exponential. 

Since  the  proportionality  constants  2//3  ,  /ir/2  ,  and  /2  or 
approximately  1.155,  1.253,  1.414  respectively  are  all  nearly  the  same, 
the  estimator  is  not  too  sensitive  to  the  possibility  that  Y  -  Xg 

and  e*  have  different  distributional  forms. 


3.  MRS.  A:  User's  Guide  and  Sample  Problem 


MRS.  A  consists  of  a  main  program  and  seven  subroutines.  The  functions 
of  these  components  are  described  In  Table  2. 

Input  instructions  are  briefly  documented  In  the  program  and  consist 
of  four  basic  card  types. 

The  first  card  or  card  image  contains  the  ancillary  statistics  for 
the  specific  problem  as  follows: 


:  Card: 

Card  Column 

Variable  Name 

Description 

1-5 

NOBS 

*  n  “  number  of  observations  (format  15;  i.e 

a  5  digit  integer,  right  justified) 

6-10 

IP 

“  p  *=  number  of  beta  parameters  (format  15) 

11-15 

ISAM 

"  K  «  number  cf  samples  for  the  mini-Mer.te 

Carlo  study  (format  15) 

16-26 

NSEED 

Ten  digit  random  number  <_  2147483647 

(format  Ill) 

28 

IWR.IT1 

■  1  if  the  main  results  of  estimation 

are  printed 

•  0  otherwise 

30 

IWRIT2 

■  1  if  the  intermediate  results  of 

estimation  are  printed 

■  0  otherwise 

32 

IWRIT3 

•  1  if  the  main  results  of  the  mini-Monte 

Carlo  study  are  printed 


0  otherwise 
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Card  Column 
34 


36 

38 


40 


42 


Variable  Name 


IWRIT4 


IWRIT5 


IWRIT6 


10PTN 


Description 

“  1  if  the  intermediate  results  of  the 
mini-Monte  Carlo  study  are  printed 
“  0  otherwise 

■  1  if  the  inputted  data  are  printed 
”  0  otherwise 

*  1  if  the  intermediate  steps  in  the 
determination  of  the  covariance  of  8 
are  printed 

*  0  otherwise 

■  1  if  the  e*'s  are  to  be  normally 


distributed 

* 

■  2  if  the  e  s  are  to  be  double  exponentially 
distributed 

■  3  if  the  e  s  are  to  be  uniformly  distributed 

IWT  ■  1  if  weights,  ,  are  to  be  assigned  to 

the  residuals 

■  0  if  residuals  are  not  weighted 


The  remaining  card  input  instructions  are  as  follows: 

Card  Number  Variable  Name  Description 

Second  card  group:  W^,  i«l,  NOBS  The  weights  assigned  to  the 

residuals  [format  (8F10.5)  ;  i.e., 
eight  ten  digit  numbers  with 
either  a  decimal  point  Included 
or  last  5  digits  are  assumed  to 
be  to  the  right  of  a  supplied 
decimal  point.] 
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Card  Number 

Variable  Name 

Description 

Third  card  group: 

i-1,  ....  NOBS 

The  observations  [format(8F10.5) ] 

Fourth  card  group: 

X^,  i-1,  ...,  NOBS 

The  matrix  of  beta  coefficients 

J-l.  IP 

read  in  by  rows  [format  (8F10.5)] 

--  The  user  may  also  supply  a  title  card  of  80  spaces  or  less  following 
the  fourth  card  group. 

The  size  of  the  problem  which  can  be  solved  Is  limited  only  by  the 
dimension  statements  in  MRS.  A.  Currently  these  restrict  the  size  to  be 
20  or  less  observations  (n  20),  10  or  less  parameters  (p  10)  and  100 
or  less  samples  in  the  mini-Monte  Carlo  study  (K  _<  100).  However,  expansion 
can  easily  be  accomplished  by  increasing  these  dimensions  in  the  dimension 
statements  as  documented  in  the  program. 

MRS.  A  is  written  in  Fortran  IV  language  and  is  compatible  with  Fortran 
G  and  H  and  WATFIV  language  compilers.  The  program  uses  double  precision 
arlthmatlc. 

MRS.  A  has  been  tested  on  several  problems  on  an  AMDAHL  470  V6  and 
should  be  compatible  with  all  IBM  computers.  MRS.  A  executes  small  problems 
such  as  n  8  5,  p  =  3,  K  ■  6  in  less  than  two  seconds.  Problems  of  sizes 
n  “  20,  p  -  2,  K  =  30  take  up  to  a  minute  of  execution  time. 

The  sample  input  and  sample  output  for  a  sample  problem  are  given  in 
Appendices  A  and  B,  respectively.  The  program  listing  is  given  in  Appendix  C. 
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Table  2 

Component 

MAIN 


INVERT 

CONST 

XTXINV 

RAND 

NORMAL 

DOUBLE 

UNIFORM 


Components  of  MRS.  A  and  Their  Functions 

Function 

Reads  data  and  generates  output. 

,  Performs  L^  estimation. 

Carries  out  the  Mini-Monte  Carlo  study. 

a  * 

Determines  a  and  the  estimated  covariance  of  B. 
Inverts  an  n  x  n  matrix. 

Constructs  the  least  squares  estimate  of  B,  Bq. 

T  —1 

Calculates  (X  X)  for  use  in  forming  Bq  . 

Generates  random  uniform  variable  with  range  0  to  1 

* 

Generates  a  vector  of  normally  distributed  e 
with  mean  0  and  variance  1  for  use  in  the  mini- 
Monte  Carlo  study. 

Generates  a  vector  of  double  exponentially 

* 

distributed  e  with  mean  0  and  variance  1  for 
use  in  the  mini-Monte  Carlo  study. 

Generates  a  vector  of  uniformly  distributed 
e  with  mean  0  and  variance  1  for  use  in  the 
mini-Monte  Carlo  study. 


.  ii 
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5  2 

20  1872530680  000 

0  0  0  1 

0 

♦•25700 

-l  .96300 

0.02400  - 

3.18000 

-3.98600 

1 .OOOOo 

1 .12000 

1.00000 

1 .95000 

1.00000 

3.02000 

1 .00000 

5.43000 

1 .00000 

6.59000 

THIS  IS 

AN  EXAMPLE 

PROBLEM  OF 

MRS.  A. 

WITH  NO  OPTIONAL  PRINTOUTS 

5  2  20  1872E39680  1  1  I  l  I  1  1  0 

4.25700  -1.98300  0.02400  -3.18000  -3.98600 

1.00000  1.12000 

1.00000  1.55000 

1.00000  3 . 02000 

1.00000  5.43000 

1.00000  6.59000 

THIS  IS  AN  EXAMPLE  PROdLEM  OF  MRS.  A.  WITH  ALU  OPTIONAL  PRINTOUTS. 
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APPENDIX  B.  SAMPLE  OUTPUTS 


MRS i  A  : 

MINIMIZES  SUM  OF  ABSOLUTE  RFSIDUALS 


-20- 


(Sample  output  with  no  optional 
pri ntouts ) 


THIS  PROGRAM  ESTIMATES  A  LINEAR  REGRESSION  BY  MINIMIZING 
THE  SUM  OF  THE  ABSOLUTE  RESIDUALS  -  LI  ESTIMATION. 

IN  ADDITION)  A  MINI-MONTE  CARLO  SIMULATION  GENERATES  AN 
ESTIMATED  COVARIANCE  MATRIX  FOR  THE  ESTIMATED 
REGRESSION  PARAMETERS. 

UNBIASED  ESTIMATES  OF  THE  REGRESSION 
PARAMETERS  ARE  OBTAINED  USING  THE  PROCEDURE 
DESCRIBED  IN  A  PAPER  BY  H.O.  HARTLEY  AND  R.L.  SIELKEN.JR. 

■TUO  LINEAR  PROGRAMMING  ALGORITHMS  FOR  UNBIASED  ESTIMATION 
OF  LINEAR  MODELS1)  1973)  JASA)  VOL.  68.  PAGES  639-41. 

THE  FOLLONING  PROCEDURE  DEVELOPED  BY  5 
D.N.  BOOK 
J.B.  BOOKER 
H.O.  HARTLEY 
R.L.  SIELKEN.  JR. 

INSTITUTE  OF  STATISTICS 
TEXAS  A  I  H  UNIVERSTIY 
COLLEGE  STATION.  TEXAS  77843 

INQUIRIES  AND  COMMENTS  SHOULD  BE  ADDRESSED  TO! 

ROBERT  L.  SIELKEN.  JR. 

THE  SUPPORT  OF  THE  OFFICE  OF  NAVAL  RESEARCH  IS  GRATEFULLY  ACKNOWLEDGED. 

THIS  IS  AN  EXAMPLE  PROBLEM  OF  MRS.  A.  WITH  NO  OPTIONAL  PRINTOUTS. 

NUMBER  OF  OBSERVATIONS  =  5 

NUMBER  OF  PARAMETERS  =  2 

THE  SAMPLE  SIZE  FOR  THE  MINI-MONTE  CARLO  STUDY  =  20 

USER  SUPPLIED  INITIAL  RANDOM  INTEGER!  NSEED  =  1872539680 

THE  AUXILIARY  LEAST  SQUARES  ESTIMATE.  BETAO.  OF  THE  REGRESSION  PARAMETER  VECTOR.  BETA 
LEAST  SQUARES  ESTIMATE  OF  BETA(  1)  =  3.227002 

LEAST  SQUARES  ESTIMATE  OF  BETA!  2)  =  -1.159747 

MRS.  A'S  ANSWER  5  THE  ESTIMATE  OF  THE  REGRESSION  PARAMETER  VECTOR  WHICH  MINIMIZES  THE  SUM  OF  THE  ABSOLUTE  RESIDUALS! 
LI  ESTIMATE  OF  BETA(  I)  =  3.416213 

LI  ESTIMATE  OF  BETA!  2)  =  -1.123249 

THE  RESIDUALS.  R(I). 1=1, NOBS 
2.098826 
3.208877 
0.000000 
0.496969 
0.000000 

THE  SUM  OF  THE  ABSOLUTE  RESIDUALS  =  5.804672 

THE  MAXIMUM  ABSOLUTE  RESIDUAL  =  3.203877 

MAIN  RESULTS  OF  THE  MINI-MONTE  CARLO  STUDY 
ESTIMATED  VALUE  OF  SIGMA  (SIGMA  HAT  4)  =  2.392353 

ESTIMATED  COVARIANCE  OF  THE  REGRESSION  PARAMETER  VECTOR  (BETA)  USING  THIS  ESTIMATE  OF  SIGMA 


4.433552 

-0.849803 


-0.849803 

0.364199 


HRS.  A  : 

MINIMIZES  SUH  OF  ABSOLUTE  RESIDUALS. 

THIS  PROGRAM  ESTIMATES  A  LINEAR  REGRESSION  BY  MINIMIZING  KamnlP  nut r.ut  w 

THE  SUM  OF  THE  ABSOLUTE  RESIDUALS  -  LI  ESTIMATION.  P  P 

IN  ADDITION.  A  MINI-MONTE  CARLO  SIMULATION  GENERATES  AN 
ESTIMATED  COVARIANCE  MATRIX  FOR  THE  ESTIMATED 
REGRESSION  PARAMETERS. 

UNBIASED  ESTIMATES  OF  THE  REGRESSION 
PARAMETERS  ARE  OBTAINED  USING  THE  PROCEDURE 
DESCRIBED  IN  A  PAPER  BY  H.O.  HARTLEY  AND  R.L.  SIELKEN.JR. 

■TWO  LINEAR  PROGRAMMING  ALGORITHMS  FOR  UNBIASED  ESTIMATION 
OF  LINEAR  MODELS* »  1973.  JASA.  VOL.  66.  PAGES  639-41. 

THE  FOLLOWING  PROCEDURE  DEVELOPED  BY  J 
D.N.  BOOK 
J.B.  BOOKER 
H.O.  HARTLEY 
R.L.  SIELKEN.  JR. 

INSTITUTE  Of  STATISTICS 

TEXAS  A  1  M  UNI VERST I Y 

COLLEGE  STATION,  TEXAS  77843 

INQUIRIES  AND  COMMENTS  SHOULD  BE  ADDRESSED  TO! 

ROBERT  L.  SIELKEN,  JR. 

THE  SUPPORT  OF  THE  OFFICE  OF  NAVAL  RESEARCH  IS  GRATEFULLY  ACKNOWLEDGED. 

THIS  IS  AN  EXAMPLE  PROBLEM  OF  MRS.  A.  WITH  ALL  OPTIONAL  PRINTOUTS. 

NUMBER  OF  OBSERVATIONS  =  5 

NUMBER  OF  PARAMETERS  =  2 

THE  SAHPLE  SIZE  FOR  THE  MINI-MONTE  CARLO  STUDY  =  20 

USER  SUPPLIED  INITIAL  RANDOM  INTEGER!  NSEED  =  1872539680 
THE  LINEAR  REGRESSION  IS  Y  =  XIBETA  F  EPSILON 
WHERE 

Y  CONTAINS  THE  OBSERVATIONS, 

BETA  IS  A  VECTOR  CONTAINING  THE  REGRESSION  PARAMETERS, 

THE  I-TH  ROW  OF  X  CONTAINS  THE  COEFFICIENTS  OF  BETA 

CORRESPONDING  TO  THE  I-TH  OBSERVATION,  AND 

EPSILCN  IS  A  RANDOM  VARIABLE  WITH  MEAN  ZERO  AND  VARIANCE 

SIGMA-SQUARED  REPRESENTING  RANDOM  VARIABILITY  FROM 

XtBETA. 

THE  Y  VECTOR 

4.25700 

-1.98300 

0.02400 

-3.18000 

-3.98600 

THE  X  MATRIX 

1.0  1.1 
1.0  2.0 
1.0  3.0 

1.0  5.4 

1.0  6.6 

SUPPLEMENTAL  INFORMATION  FROM  THE  LINEAR  PROGRAMMING  PROBLEM  DETERMINATION  OF  THE  ESTIMATE  OF 
THE  REGRESSION  PARAMETER  VECTOR,  BETA. 

THE  LINEAR  PROGRAMMING  PROBLEM  AS  IT  WAS  CREATED 

THE  OBJECTIVE  FUNCTION  COEFFICIENTS 
C(  1)  =  0.00000 

C<  2)  =  0.00000 

C<  3)  =  0.00000 

C<  4)  =  0.00000 

C<  5)  =  -1.00000 

C<  6)  =  -1.00000 

C<  7)  =  -1,00000 

C(  8)  *  -1.00000 

C<  9)  =  -1.00000 

C<  10)  =  0.00000 

C<  11)  =  0.00000 

C(  12)  =  0.00000 

C<  13)  =  0.00000 

C<  14)  =  0.00000 

C<  15)  =  0.00000 

C<  16)  =  0.00000 

C(  17)  =  0.00000 

C(  18)  =  0.00000 

C(  19)  =  0.00000 

THE  CONSTRAINT  MATRIX  A 

-1.0  -1.1  1.0  1.1  -1.0  0.0  0.0  0.0  0.0  1.0  0.0  0.0  0.0 


(Sample  output  with  optional  printouts) 


THE  VARIABLES  INITIALLY  IN  THE  BASIS 
10,11,12,13,14,15,16,17,18,19, 

THE  PROBLEM  P2  HAS  BEEN  SELECTED 
THE  INVERSE  OF  XTX 

0.810298  -0.168497 

-0.168497  0.046521 

VALUE  OF  BETAO  TO  COMPUTE  RHS 
3.227002  -1.159747 

THE  RHS,  YMXB,  FOR  PI  OR  P2 
2.328914 
-2.948496 
0.299433 
-0.109578 
0.429728 

THE  INITIAL  VALUES  OF  THE  BASIC  VARIABLES 
XB(  1)  =  0.23289D  01 

XB(  2)  =  -0.29485D  01 

XB<  3)  =  0.2994 3D  00 

XB(  4)  =  -0.10958D  00 

XB(  5)  =  0.42973D  00 


-0.23289D  01 
0.29485D  01 
-0.29943D  00 
0.10958D  00 
-0.42973D  00 
THE  INITIAL  VALUE  OF  THE  OBJECTIVE  FUNCTION  = 
THE  REDUCED  COSTS 


XB( 

XB( 

XB( 

XB( 

XB( 


6)  = 

7)  = 

8)  = 
9)  = 

10)  = 


O.OOOOOD  00 


THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 


1- TH  REDUCED  COST 

2- TH  REDUCED  COST 

3- TH  REDUCED  COST 

4- TH  REDUCED  COST 

5- TH  REDUCED  COST 

6- TH  REDUCED  COST 

7- TH  REDUCED  COST 

8- TH  REDUCED  COST 

9- TH  REDUCED  COST 

10- TH  REDUCED  COST 

1 1- TH  REDUCED  COST 

12- TH  REDUCED  COST 

13- TH  REDUCED  COST 

14- TH  REDUCED  COST 

15-  ,H  REDUCED  COST 

1 6- TH  REDUCED  COST 

17- TH  REDUCED  COST 

18- TH  REDUCED  COST 

19- TH  REDUCED  COST 


-O.OOOOOD  00 
-O.OOOOOD  00 
-O.OOOOOD  00 
-O.OOOOOD  00 
-0.10000D  01 
-0.10000D  01 
-0.10000D  01 
-0.10000D  01 
-0.10000D  01 
-O.OOOOOD  00 
-O.OOOOOD  00 
-O.OOOOOD  00 
-O.OOOOOD  00 
-O.OOOOOD  00 
-O.OOOOOD  00 
-O.OOOOOD  00 
-O.OOOOOD  00 
-O.OOOOOD  00 
-O.OOOOOD  00 


THE  ... 

TENTATIVELY  THE 
THE  YRJ'S 

(  YR(J)  =  0  IF  THE  J-TH  VARIABLE  IS  IN  THE  BASIS. 


2-TH  BASIC  VARIABLE  IS  LEAVING  THE  BASIS 


)N  BE  ZERO 

OTHERWISE  TOO.) 

YR! 

1)  = 

-0.10000D  01 

YR( 

2)  = 

-0.19500D  01 

YR( 

3)  = 

0.10000D  01 

YR( 

4)  = 

0.19500D  01 

YR( 

5)  = 

O.OOOOOD  00 

YR( 

6)  = 

-0.10000D  01 

YR< 

7)  = 

O.OOOOOD  00 

YR( 

8)  = 

O.OOOOOD  00 

YR( 

9)  = 

O.OOOOOD  00 

YR< 

10)  = 

O.OOOOOD  00 

YR< 

11)  = 

O.OOOOOD  00 

YR( 

12)  = 

O.OOOOOD  00 

YR (  13)  =  O.OOOOOD  00 

YRt  14)  =  O.OOOOOD  00 

YRt  15)  =  O.OOOOOD  00 

YRt  16)  =  O.OOOOOD  00 

YR(  17)  =  O.OOOOOD  00 

YRt  18)  =  O.OOOOOD  00 

YRt  19)  =  O.OOOOOD  00 

THE  RATIO  OF  THE  NET  PRICES  TO  THE  NEGATIVE  YRJ'S 

THE  RATIO  Of  THE  1-TH  NET  PRICE  TO  YRt  1) 

THE  RATIO  OF  THE  2-TH  NET  PRICE  TO  YRt  2) 

THE  RATIO  OF  THE  6-TH  NET  PRICE  TO  YRt  0) 

THE  2-TH  VARIABLE  IS  LEAVING  THE  BASIS. 

THE  2-TH  VARIABLE  IS  ENTERING  THE  BASIS. 

THE  BASIC  VARIABLES  ARE  NON 

10>  2»  12>  13i  14>  15)  16)  17>  18) 

THE  VALUES  OF  THE  BASIC  VARIABLES  ARE  NON 
XBt  1)  =  0.40224D  01 

XBt  2)  =  0.15120D  01 

XBt  3)  =  0.48658D  01 

XBt  4)  =  0.8J009D  01 

XBt  5)  =  0.10294D  02 

XBt  6)  =  -0.40224D  01 

XBt  7)  =  0.22204D-15 

XBt  8)  =  -0.48658D  01 

XBt  9)  =  -0.81009D  01 

XBt  10)  =  -0.10394D  02 

THE  CURRENT  VALUE  OF  THE  OBJECTIVE  FUNCTION  = 

THE  REDUCED  COSTS 

THE  1-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  2-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  3-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  4-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  5-TH  REDUCED  COST  =  -0.10000D  01 

THE  6-TH  REDUCED  COST  =  -0.10000D  01 

THE  7-TH  REDUCED  COST  =  -0.10000D  01 

THE  8-TH  REDUCED  COST  =  -0.10000D  01 

THE  9-TH  REDUCED  COST  =  -0.10000B  01 

THE  10-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  11-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  12-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  13-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  14-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  15-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  16-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  17-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  18-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  19-TH  REDUCED  COST  =  -O.OOOOOD  00 


-O.OOOOOD  00 
-O.OOOOOD  00 
-0.10000D  01 


18)  19i 


XBt  10)  = 


-O.OOOOOD  00 


-O.OOOOOD  00 
-O.OOOOOD  00 
-O.OOOOOD  00 
-O.OOOOOD  00 
-0.10000D  01 
-0.10000D  01 
-0.10000D  01 
-0.10000D  01 
-0.10000D  01 
-O.OOOOOD  00 
-O.OOOOOD  00 
-O.OOOOOD  00 
-O.OOOOOD  00 
-O.OOOOOD  00 
-O.OOOOOD  00 
-O.OOOOOD  00 
-O.OOOOOD  00 
-O.OOOOOD  00 
-O.OOOOOD  00 


TENTATIVELY  THE  10-TH  BASIC  VARIABLE  IS  LEAVING  THE  BASIS 
THE  YRJ'S 

(  YRt J)  =  0  IF  THE  J-TH  VARIABLE  IS  IN  THE  BASIS. 


IT  CAN  BE  ZERO  OTHERWISE  TOO.) 

YRt  1) 

= 

-0.23795D  01 

YRt  2) 

r 

O.OOOOOD  00 

YRt  3) 

= 

0.23795D  01 

YRt  4) 

= 

-0.44409D-15 

YRt  5) 

= 

O.OOOOOD  00 

YRt  6) 

= 

-0.33795D  01 

YRt  7) 

s 

O.OOOOOD  00 

YRt  8) 

2 

O.OOOOOD  00 

YRt  9) 

= 

-0.10000D  01 

YRt  10) 

= 

O.OOOOOD  00 

YRt  11) 

= 

0.33795D  01 

YRt  12) 

= 

O.OOOOOD  00 

YRt  13) 

r 

O.OOOOOD  00 

YRt  14) 

s 

O.OOOOOD  00 

YRt  15) 

r 

O.OOOOOD  00 

YRt  16) 

s 

O.OOOOOD  00 

YRt  17) 

s 

O.OOOOOD  00 

YRt  18) 

s 

O.OOOOOD  00 

YRt  19) 

s 

O.OOOOOD  00 

TO  RATIO  OF  TO  NET  PRICES  TO  THE  NEGATIVE  YRJ'S 

THE  RATIO  OF  THE  1-TH  NET  PRICE  TO  YRt 

THE  RATIO  OF  THE  6-TH  NET  PRICE  TO  YRt 

THE  RATIO  OF  THE  9-TH  NET  PRICE  TO  YRt 

TO  10-TH  VARIABLE  IS  LEAVING  THE  BASIS. 

THE  1-TH  VARIABLE  IS  ENTERING  THE  BASIS. 

THE  BASIC  VARIABLES  ARE  NOV 

10)  2)  12*  13)  14»  15)  16i  17) 

TO  VALUES  Of  THE  BASIC  VARIABLES  ARE  NOW 
XBt  1)  *  0.58817D  01 

XBt  2)  *  -0.72807B  00 


-O.OOOOOD  00 
-0.29590D  00 
-0.10000D  01 
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XB(  3)  *  0.246B9D  01 

XB(  4)  =  0.30525D  00 

XBI  5)  =  0 . 6661 3D-15 

XBI  6)  =  -0.58817D  01 

XB(  7)  =  0. 100800-15 

X8(  8)  =  -0.24689D  01 

XBI  9)  =  -0.30525D  00 

XBI  10)  =  0.43682D  01 

THE  CURRENT  VALUE  OF  THE  OBJECTIVE  FUNCTION  =  -O.OOOOOD  00 
THE  REDUCED  COSTS 

THE  1-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  2-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  3-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  4-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  5-TH  REDUCED  COST  =  -0.10000D  01 

THE  6-TH  REDUCED  COST  =  -0.10000D  01 

THE  7-TH  REDUCED  COST  =  -0.10000P  01 

THE  8-TH  REDUCED  COST  =  -0.10000D  01 

THE  9-TH  REDUCED  COST  =  -0.10000D  01 

THE  10-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  11-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  12- TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  13-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  14-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  15-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  16-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  17-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  18-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  19-TH  REDUCED  COST  =  -O.OOOOOD  00 

TENTATIVELY  THE  6-TH  BASIC  VARIABLE  IS  LEAVING  THE  BASIS 
THE  YRJ'S 

(  '  YR( J)  =  0  IF  THE  J-TH  VARIABLE  IS  IN  THE  BASIS. 

IT  CAN  BE  ZERO  OTHERWISE  TOO,) 

YRI  1)  =  O.OOOOOD  00 

YR(  2)  =  O.OOOOOD  00 

YR(  3)  =  -0.18041D-15 

YR(  4)  =  -0.44409D-15 

YR(  5)  =  -0.10000D  01 

YR(  6)  =  -0.11789D  01 

YR(  7)  =  O.OOOOOD  00 

YRI  8)  =  O.OOOOOD  00 

YR(  9)  =  -0.17888D  00 

YR(  10)  *  O.OOOOOD  00 

YR(  11)  =  0.11 7890  01 

YR(  12)  =  O.OOOOOD  00 

YR(  13)  =  O.OOOOOD  00 

YR{  14)  =  O.OOOOOD  00 

YR(  15)  =  O.OOOOOD  00 

YR(  16)  =  O.OOOOOD  00 

YR(  17)  =  O.OOOOOD  00 

YR(  18)  =  O.OOOOOD  00 

YR(  19)  =  0.17888D  00 

THE  RATIO  OF  THE  NET  PRICES  TO  THE  NEGATIVE  YRJ'S 

THE  RATIO  OF  THE  5-TH  NET  PRICE  TO  TR(  5)  =  -0.10000D  01 

THE  RATIO  OF  THE  6-TH  NET  PRICE  TO  YR(  6)  =  -0.B4826D  00 

THE  RATIO  OF  TIC  9-TH  NET  PRICE  TO  YR(  9)  =  -0.55904D  01 

THE  6-TH  VARIABLE  IS  LEAVING  THE  BASIS. 

THE  6-TH  VARIABLE  IS  ENTERING  THE  BASIS. 

THE  BASIC  VARIABLES  ARE  NON 


10)  2*  12*  13i  14j  6»  16*  17)  18)  1) 

THE  VALUES  OF  THE  BASIC  VARIABLES  ARE  NON 


XB( 

1)  = 

0.22204D-15 

XB( 

2)  = 

0.34720D  00 

XB< 

3)  = 

-0. 136980  01 

XBI 

4)  = 

-0.942060  00 

XBI 

5)  = 

-0.154950-14 

XB( 

6)  = 

0.49892D  01 

XBI 

7)  = 

0.99785D  01 

XBI 

8)  = 

0.136980  01 

XBI 

9)  = 

0.94206D  00 

XBI 

10)  = 

-0.271780  01 

THE  CURRENT  VALUE  OF  THE  OBJECTIVE  FUNCTION  *  -0.49892D  01 

THE  REDUCED  COSTS 

THE  i-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  2-TH  REDUCED  COST  =  -O.OOOOOD  00 


THE  3-TH  REDUCED  COST  * 
THE  4-TH  REDUCED  COST  = 
THE  5-TH  REDUCED  COST  = 
THE  6-TH  REDUCED  COST  = 


0.I3B78D-15 
0.27756D-15 
-0.151740  00 
-O.OOOOOD  00 


THE  7-TH  REDUCED  COST  *  -O.IOOOOD  01 


THE  8-TH  REDUCED  COST  * 


-O.IOOOOD  01 
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THE  9-TH  REDUCED  COST 
THE  10-TH  REDUCED  COST 
THE  11 -TH  REDUCED  COST 
THE  12-TH  REDUCED  COST 
THE  13-TH  REDUCED  COST 
THE  H-TH  REDUCED  COST 
THE  15-TH  REDUCED  COST 
THE  14-TH  REDUCED  COST 
THE  17-TH  REDUCED  COST 
THE  18-TH  REDUCED  COST 
THE  19-TH  REDUCED  COST 


-0.84824D  00 
-O.OOOOOD  00 
-0.10000D  01 
-O.OOOOOD  00 
-O.OOOOOD  00 
-O.OOOOOD  00 
-0.84824D  00 
-O.OOOOOD  00 
-O.OOOOOD  00 
-O.OOOOOD  00 
-0.15174D  00 


THE  19-TH  REDUCED  COST  =  -0.15174D  00 

TENTATIVELY  THE  10-TH  BASIC  VARIABLE  IS  LEAVING  THE  BASIS 
THE  YRJ'S 

(  YR(J)  =  0  IF  THE  J-TH  VARIABLE  IS  IN  THE  BASIS. 


IT  CAN  BE  ZERO  OTHERWISE  TOO.) 

YR( 

1) 

S 

O.OOOOOD  00 

YR( 

2) 

r 

O.OOOOOD  00 

YRC 

3) 

s 

-0.10000D  01 

YR( 

4) 

= 

O.OOOOOD  00 

YR( 

5) 

= 

-0.1 204 8D  01 

YR( 

4) 

= 

O.OOOOOD  00 

YR( 

7) 

r 

O.OOOOOD  00 

YR( 

8) 

= 

O.OOOOOD  00 

YR( 

9) 

= 

0.20475D  00 

YR( 

10) 

s 

O.OOOOOD  00 

YR( 

11) 

s 

O.OOOOOD  00 

YR( 

12) 

_ 

O.OOOOOD  00 

YR( 

13) 

= 

O.OOOOOD  00 

YRI 

14) 

= 

O.OOOOOD  00 

YR( 

15) 

= 

0.12048D  01 

YRI 

14) 

= 

O.OOOOOD  00 

YRI 

17) 

= 

O.OOOOOD  00 

YRI 

18) 

= 

O.OOOOOD  00 

YRI 

19) 

s 

-0.20475D  00 

THE  RATIO  OF  THE  NET  PRICES  TO  THE  NEGATIVE  YRJ'S 

THE  RATIO  OF  THE  3-TH  NET  PRICE  TO  YRI  3)  = 

THE  RATIO  OF  THE  5-TH  NET  PRICE  TO  YR<  5)  = 

THE  RATIO  OF  THE  19-TH  NET  PRICE  TO  YRI  19)  = 

THE  10-TH  VARIABLE  IS  LEAVING  THE  BASIS. 

THE  3-TH  VARIABLE  IS  ENTERING  THE  BASIS. 

THE  BASIC  VARIABLES  ARE  NON 

10.  2.  12.  13.  14.  4>  16,  17,  18, 

THE  VALUES  OF  THE  BASIC  VARIABLES  ARE  NON 
XB<  1)  =  0.14441D-15 

XB(  2)  =  0.34720D  00 

XBI  3)  =  -0.1349BD  01 

XB(  4)  =  -0.94206D  00 

XB(  5)  =  -0.531 18D-15 

XB(  6)  =  0.49892D  01 

XB(  7)  =  0.99785D  01 

XB(  8)  =  0.13698D  01 

XB<  9)  =  0.94204D  00 

XB(  10)  =  0.27178D  01 

THE  CURRENT  VALUE  OF  THE  OBJECTIVE  FUNCTION  =  -0. 

THE  REDUCED  COSTS 


0.13878D-15 
-0.12595D  00 
-0.74107D  00 


XB(  10)  = 


-0.49892D  01 


THE  1-TH 
THE  2-TH 
THE  3-TH 
THE  4-TH 
THE  5-TH 
THE  4-TH 
THE  7-TH 
THE  8-TH 
THE  9-TH 
THE  10-TH 
THE  11-TH 
THE  12-TH 
THE  13-TH 
THE  14-TH 
THE  15-TH 
THE  14-TH 
THE  17-TH 
THE  18-TH 
THE  19-TH 
TENTATIVELY  THE 
THE  YRJ'S 
(  YR(J)  =  0  IF 


1- TH  REDUCED  COST  = 

2- TH  REDUCED  COST  = 

3- TH  REDUCED  COST  = 

4- TH  REDUCED  COST  = 

5- TH  REDUCED  COST  - 
4-TH  REDUCED  COST  = 

7- TH  REDUCED  COST  = 

8- TH  REDUCED  COST  = 

9- TH  REDUCED  COST  = 

10- TH  REDUCED  COST  = 

11- TH  REDUCED  COST  = 

12- TH  REDUCED  COST  = 

13- TH  REDUCED  COST  = 

14- TH  REDUCED  COST  = 

15- TH  REDUCED  COST  = 
14-TH  REDUCED  COST  = 

17- TH  REDUCED  COST  = 

18- TH  REDUCED  COST  = 

19- TH  REDUCED  COST  = 


-O.OOOOOD  00 
-O.OOOOOD  00 
-O.OOOOOD  00 
0.24368D-15 
-0.15174D  00 
-O.OOOOOD  00 
-0.10000D  01 
-0.10000D  01 
-0.84824D  00 
-O.OOOOOD  00 
-0.10000D  01 
-O.OOOOOD  00 
-O.OOOOOD  00 
-O.OOOOOD  00 
-0.84824D  00 
-O.OOOOOD  00 
-O.OOOOOD  00 
-O.OOOOOD  00 
-0.15174D  00 


3-TH  BASIC  VARIABLE  IS  LEAVING  THE  BASIS 


0  IF  THE  J-TH  VARIABLE  IS  IN  THE  BASIS. 


IT  CAN  BE  ZERO  OTHERWISE  TOO.) 
YR(  1)  =  -0.13878B-14 

YRI  2)  =  O.OOOOOD  00 

YR<  3)  =  O.OOOOOD  00 


•t-  . 
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YR(  4)  =  O.OOOOOD  00 

YR(  5)  =  -0.A5265D  00 

YRI  6)  =  O.OOOOOD  00 

YR(  7)  =  -0.10000D  01 

YR<  8)  =  O.OOOOOD  00 

YR(  9)  =  -0.34735D  00 

YR(  10)  =  O.OOOOOD  00 

YRI  11)  =  O.OOOOOD  00 

YR(  12)  =  O.OOOOOD  00 

YR(  13)  =  O.OOOOOD  00 

YR(  14)  =  O.OOOOOD  00 

YR(  15)  =  0.65265D  00 

YR(  16)  =  O.OOOOOD  00 

YRI  17)  =  O.OOOOOD  00 

YR(  18)  =  O.OOOOOD  00 

YR(  19)  =  0.34735D  00 

THE  RATIO  OF  THE  NET  PRICES  TO  THE  NEGATIVE  YRJ'S 

THE  RATIO  OF  THE  5-TH  NET  PRICE  TO  YR(  5)  =  -0.23249D  00 

THE  RATIO  OF  THE  7-TH  NET  PRICE  TO  YR(  7)  =  -0.10000D  01 

THE  RATIO  OF  THE  9-TH  NET  PRICE  TO  YR(  9)  =  -0.24421D  01 

THE  3-TH  VARIABLE  IS  LEAVING  THE  BASIS. 

THE  5-TH  VARIABLE  IS  ENTERING  THE  BASIS. 

THE  BASIC  VARIABLES  ARE  NOW 

10.  2i  5,  13.  H»  6,  16>  17,  18,  3, 

THE  VALUES  OF  THE  BASIC  VARIABLES  ARE  NOW 
XB<  1)  =  0.41977D  01 

XB<  2)  =  -0.3&497D-01 

XB(  3)  =  0.20988D  01 

XB<  4)  =  -0.49697D  00 

XB<  5)  =  -0.68800D-15 

XB(  6)  =  0.32089D  01 

XBI  7)  =  0.64178D  01 

XBI  8)  =  0.22204D-15 

XBI  9)  =  0.49697D  00 

XBI  10)  =  0.18921D  00 

THE  CURRENT  VALUE  OF  THE  OBJECTIVE  FUNCTION  =  -0.53077D  01 

THE  REDUCED  COSTS 

THE  1-TH  REDUCED  COST  =  -0.20817B-15 

THE  2-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  3-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  4-TH  REDUCED  COST  =  0.22204D-15 

THE  5-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  6-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  7-TH  REDUCED  COST  =  -0.76751D  00 

THE  8-TH  REDUCED  COST  *  -0.10000D  01 

THE  9-TH  REDUCED  COST  =  -0.76751D  00 

THE  10-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  11-TH  REDUCED  COST  =  -Q.10000D  01 

THE  12-TH  REDUCED  COST  =  -0.23249D  00 

THE  13-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  14-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  15-TH  REDUCED  COST  =  -0.10000D  01 

TIE  14-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  17-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  18-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  19-TH  REDUCED  COST  =  -0.23249D  00 

TENTATIVELY  THE  4-TH  BASIC  VARIABLE  IS  LEAVING  THE  BASIS 
THE  YRJ'S 

I  YRIJ)  =  0  IF  THE  J-TH  VARIABLE  IS  IN  THE  BASIS. 


IT  CAN  BE  ZERO  OTHERWISE  TOO.) 

YRI 

1) 

= 

-0.13878D-16 

YRI 

2) 

= 

O.OOOOOD  00 

YRI 

3) 

s 

O.OOOOOD  00 

YRI 

4) 

= 

-0.66613D-15 

YRI 

5) 

= 

O.OOOOOD  00 

YRI 

6) 

s 

O.OOOOOD  00 

YRI 

7) 

= 

0.32493D  00 

YRI 

8) 

= 

-0.10000D  01 

YRI 

9) 

z 

-0.67507D  00 

YRI 

10) 

s 

O.OOOOOD  00 

YRI 

11) 

z 

O.OOOOOD  00 

YRI 

12) 

z 

-0.32493D  00 

YRI 

13) 

z 

O.OOOOOD  00 

YRI  14) 

z 

O.OOOOOD  00 

YRI  15) 

z 

O.OOOOOD  00 

YRI  14) 

z 

O.OOOOOD  00 

YRI 

17) 

z 

O.OOOOOD  00 

YRI  18) 

z 

O.OOOOOD  00 

YRI 

19) 

z 

0.47507D  00 

TIC  RATIO  OF  THE  NET  PRICES  TO  THE  NEGATIVE  YRJ'S 

THE  RATIO  OF  THE  8-TH  NET  PRICE  TO  YR(  8)  =  -0.10000D  01 


THE  RATIO  Of  TK  9-TH  MET  PRICE  TO  YR<  9)  *  -0.11369D  01 

THE  RATIO  OF  THE  12-TH  NET  PRICE  TO  YR(  12)  =  -0.71552D  00 

THE  4-TH  VARIABLE  IS  LEAVING  THE  BASIS. 

THE  12-TH  VARIABLE  IS  ENTERING  THE  BASIS. 

THE  BASIC  VARIABLES  ARE  NOW 

lOi  2t  5.  12>  14>  6)  l&i  17)  18>  3> 

THE  VALUES  OF  THE  BASIC  VARIABLES  ARE  NOU 
XB<  1)  =  0.88846D  01 

XB(  2)  =  -0.46492D  00 

XB<  3)  =  0.44423D  01 

XB(  4)  =  0.15295D  01 

XB(  5)  =  -0.86309D-15 

XB(  6)  =  0.12210D  01 

XB<  7)  =  0.24420D  01 

XB<  8)  =  -0.15295D  01 

XB<  9)  =  0.555110-16 

XB(  10)  =  -0.263410  01 

THE  CURRENT  VALUE  OF  THE  OBJECTIVE  FUNCTION  =  -0.566330  01 

THE  REDUCED  COSTS 

THE  1-TH  REDUCED  COST  =  -0.124900-15 

THE  2-TH  REDUCED  COST  =  -0.000000  00 

THE  3-TH  REDUCED  COST  =  -0.000000  00 

THE  4-TH  REDUCED  COST  =  0.666130-15 

THE  5-TH  REDUCED  COST  =  -O.OOOOOB  00 

THE  6-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  7-TH  REDUCED  COST  =  -0.10000D  01 

THE  8-TH  REDUCED  COST  =  -0.28448D  00 

THE  9-TH  REDUCED  COST  =  -0.28448D  00 

THE  10-TH  REDUCED  COST  =  -0.000000  00 

THE  1 1-TH  REDUCED  COST  -  -0.10000D  01 

THE  12-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  13-TH  REDUCED  COST  =  -0.71552D  00 

THE  14-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  15-TH  REDUCED  COST  =  -0.10000D  01 

THE  16-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  17-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  18-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  19-TH  REDUCED  COST  =  -0.71552D  00 

TENTATIVELY  THE  10-TH  BASIC  VARIABLE  IS  LEAVING  THE  BASIS 
THE  YRJ'S 

(  YR(J)  =  0  IF  THE  J-TH  VARIABLE  IS  IN  THE  BASIS. 


IT  CAN  BE  ZERO  OTHERWISE  TOO.) 

YR( 

1) 

= 

-0.10000D  01 

YR( 

2) 

= 

O.OOOOOD  00 

YR( 

3) 

= 

O.OOOOOD  00 

YR( 

4) 

= 

-0.35527D-14 

YR( 

5) 

= 

O.OOOOOD  00 

YR( 

6) 

= 

O.OOOOOD  00 

YR( 

7) 

S 

O.OOOOOD  00 

YR( 

8) 

= 

-0.56B10D  01 

YR( 

9) 

= 

-0.46810D  01 

YR( 

10) 

= 

O.OOOOOD  00 

YR( 

11) 

s 

O.OOOOOD  00 

YR( 

12) 

= 

O.OOOOOD  00 

YR( 

13) 

= 

0.56810D  01 

YR< 

14) 

= 

O.OOOOOD  00 

YR( 

15) 

= 

O.OOOOOD  00 

YR< 

16) 

= 

O.OOOOOD  00 

YR( 

17) 

= 

O.OOOOOD  00 

YR( 

18) 

= 

O.OOOOOD  00 

YR( 

19) 

= 

0.46810D  01 

THE  RATIO  OF  THE  NET  PRICES  TO  THE  NEGATIVE  YRJ'S 

THE  RATIO  OF  THE  1-TH  NET  PRICE  TO  YR(  1)  =  -0.12490D-15 

THE  RATIO  OF  THE  8-TH  NET  PRICE  TO  YR(  8)  =  -0.50076D-01 

THE  RATIO  OF  THE  9-TH  NET  PRICE  TO  YR(  9)  =  -0.60773D-01 

THE  10-TH  VARIABLE  IS  LEAVING  THE  BASIS, 

THE  1-TH  VARIABLE  IS  ENTERING  THE  BASIS. 

THE  BASIC  VARIABLES  ARE  NOU 

10)  2>  5i  12)  14*  6>  16)  17)  18>  1> 

THE  VALUES  OF  THE  BASIC  VARIABLES  ARE  NOU 
XB(  1)  =  0.88846D  01 

XB(  2)  *  — 0.46492D  00 

XB<  3)  =  0.44423D  01 

XB(  4)  =  0.15295D  01 

XB<  5)  =  -0.86309D-15 

XB(  6)  =  0.12210D  01 

XB(  7)  =  0.24420D  01 

XB(  8)  =  -0.15295D  01 

XB<  9)  =  0.18956D-16 

XB(  10)  =  0.26341D  01 

THE  CURRENT  VALUE  OF  THE  OBJECTIVE  FUNCTION  =  -0.56633D  01 
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THE  REDUCED  COSTS 

THE  1-TH  REDUCED  COST 
THE  2-TH  REDUCED  COST 
THE  3-TH  REDUCED  COST 
THE  4-TH  REDUCED  COST 
THE  5-TH  REDUCED  COST 
THE  4-TH  REDUCED  COST 
THE  7-TH  REDUCED  COST 
THE  8-TH  REDUCED  COST 
THE  9-TH  REDUCED  COST 
THE  10-TH  REDUCED  COST 
THE  11-TH  REDUCED  COST 
THE  12-TH  REDUCED  COST 
THE  13-TH  REDUCED  COST 
THE  14-TH  REDUCED  COST 
THE  15-TH  REDUCED  COST 
THE  14-TH  REDUCED  COST 
THE  17-TH  REDUCED  COST 
THE  18-TH  REDUCED  COST 
19-TH  REDUCED  COST 


-O.OOOOOD  00 
-O.OOOOOD  00 
0.41433D-14 
0.44413D-15 
-O.OOOOOD  00 
-O.OOOOOD  00 
-0.10000D  01 
-0.28448D  00 
-0.2B448D  00 
-O.OOOOOD  00 
-0.10000D  01 
-O.OOOOOD  00 
-0.71552D  00 
-O.OOOOOD  00 
-0.10000D  01 
-O.OOOOOD  00 
-O.OOOOOD  00 
-O.OOOOOD  00 
-0.71552D  00 


8-TH  BASIC  VARIABLE  IS  LEAVING  THE  BASIS 


-0.10000D  01 
-0.92437D-01 
-0.13493D  00 


THE  ..  ... 

TENTATIVELY  THE 

THE  YRJ'S 

<  WJ)  =  0  If  THE  J-TH  VARIABLE  IS  IN  THE  BASIS. 

IT  CAN  BE  ZERO  OTHERWISE  TOO.) 

YR<  1)  =  O.OOOOOD  00 

YR(  2)  =  O.OOOOOD  00 

YR<  3)  =  O.OOOOOD  00 

YR(  4)  =  -0.15543D-14 

YR(  5)  =  O.OOOOOD  00 

YR(  4)  =  O.OOOOOD  00 

YR(  7)  =  -0.10000D  01 

YR<  8)  =  -0.30774D  01 

YR<  9)  =  -0.20774D  01 

YR(  10)  =  O.OOOOOD  00 

YR<  11)  =  O.OOOOOD  00 

YR(  12)  =  O.OOOOOD  00 

YR(  13)  =  0.30774D  01 

YR<  14)  =  O.OOOOOD  00 

YR(  15)  =  O.OOOOOD  00 

YR<  14)  =  O.OOOOOD  00 

YR(  17)  =  O.OOOOOD  00 

YR(  18)  =  O.OOOOOD  00 

YR(  19)  =  0.20774D  01 

THE  RATIO  OF  THE  NET  PRICES  TO  THE  NEGATIVE  YRJ'S 

THE  RATIO  OF  THE  7-TH  NET  PRICE  TO  YR(  7)  = 

THE  RATIO  OF  THE  8-TH  NET  PRICE  TO  YR  8)  = 

THE  RATIO  OF  THE  9-TH  NET  PRICE  TO  YR(  9)  = 

THE  8-TH  VARIABLE  IS  LEAVING  THE  BASIS. 

THE  8-TH  VARIABLE  IS  ENTERING  THE  BASIS. 

THE  BASIC  VARIABLES  ARE  NOW 

10.  2i  5i  12i  14»  4>  14>  8>  1B»  1> 

THE  VALUES  OF  THE  BASIC  VARIABLES  ARE  NOW 
0.41977D  01 
-0.34497D-01 
0.20988D  01 
0.44409D-15 
-0.48800D-15 
0.32089D  01 
0.44178D  01 
0.49497D  00 
0.99394D  00 
.  -0.18921D  00 

« SSS  Kf  " ,HE rm,m  =  « 

THE  1-TH  REDUCED  COST  =  -O.OOOOOD  00 

2- TH  REDUCED  COST  *  -O.OOOOOD  00 

3- TH  REDUCED  COST  =  -0.12490D-15 

4- TH  REDUCED  COST  =  0.44413D-15 

5- TH  REDUCED  COST  =  -O.OOOOOD  00 

4-TH  REDUCED  COST  *  -O.OOOOOD  00 

7- TH  REDUCED  COST  *  -0.90754D  00 

8- TH  REDUCED  COST  =  -O.OOOOOD  00 

9- TH  REDUCED  COST  =  -0.92437D-01 

10- TH  REDUCED  COST  =  -O.OOOOOD  00 

11- TH  REDUCED  COST  =  -0.10000D  01 

12- TH  REDUCED  COST  *  -O.OOOOOD  00 

13- TH  REDUCED  COST  *  -0.10000D  01 

14- TH  REDUCED  COST  *  -O.OOOOOD  00 

15- TH  REDUCED  COST  =  -0.10000D  01 

14-TH  REDUCED  COST  *  -O.OOOOOD  00 

17-TH  REDUCED  COST  »  -0.92437D-01 


XB( 

1)  = 

XB< 

2)  = 

XB( 

3)  = 

X8( 

4)  = 

XB( 

5)  = 

XB( 

4)  = 

XB( 

7)  = 

XB( 

8)  = 

XB( 

9)  = 

XB( 

10)  = 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 


-29- 


TH€  18-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  19- TH  REDUCED  COST  *  -0.90756D  00 

TENTATIVELY  THE  10-TH  BASIC  VARIABLE  IS  LEAVING  THE  BASIS 
THE  YRJ'S 

(  YR ( J)  -  0  IF  THE  J-TH  VARIABLE  IS  IH  THE  BASIS. 


IT  CAN  BE  ZERO  OTHERWISE  TOO.) 

YR( 

1) 

= 

O.OOOOOD  00 

YRf 

2) 

5 

O.OOOOOD  00 

YR( 

3) 

= 

-0.10000D  01 

YRf 

4) 

= 

0.46629D-14 

YRf 

5) 

= 

O.OOOOOD  00 

YRf 

6) 

= 

O.OOOOOD  00 

YRf 

7) 

s 

-0.18459D  01 

YRf 

8) 

= 

O.OOOOOD  00 

YRf 

9) 

= 

0.84594B  00 

YRf 

10) 

s 

O.OOOOOD  00 

YRf 

11) 

s 

O.OOOOOD  00 

YRf 

12) 

= 

O.OOOOOD  00 

YRf 

13) 

O.OOOOOD  00 

YRf 

14) 

= 

O.OOOOOD  00 

YRf 

15) 

= 

O.OOOOOD  00 

YRf 

16) 

s 

O.OOOOOD  00 

YRf 

17) 

= 

0.18459D  01 

YRf 

18) 

= 

O.OOOOOD  00 

YRf 

19) 

r 

-0.84594D  00 

THE  RATIO  OF 

THE  NET  PRICES  TO 

THE  RATIO  OF  THE  3-TH  NET  PRICE  TO  YR( 

THE  RATIO  OF  THE  7-TH  NET  PRICE  TO  YR( 

THE  RATIO  OF  THE  19-TH  NET  PRICE  TO  YR< 

THE  10-TH  VARIABLE  IS  LEAVING  THE  BASIS. 

THE  3-TH  VARIABLE  IS  ENTERING  THE  BASIS. 

THE  BASIC  VARIABLES  ARE  NON 

10.  2.  5.  12.  14.  6.  16.  8. 

THE  VALUES  OF  THE  BASIC  VARIABLES  ARE  NON 


iJ'S 
3)  = 
7)  = 
19)  = 


-0.12490D-15 
-0.49165D  00 
-0.10728D  01 


18.  3. 


XBf 

2)  = 

-0.36497D-01 

XBf 

3)  = 

0.20988D  01 

XBf 

4)  = 

0.402081-15 

XBf 

5)  = 

-0.69062D-15 

XBf 

6)  = 

0.32089D  01 

XBf 

7)  = 

0.64178D  01 

XBf 

8)  = 

0.49697D  00 

XBf 

9)  = 

0.99394D  DO 

XBf 

10)  = 

0.18921D  00 

THE  CURRENT  VALUE  OF  THE  OBJECTIVE  FUNCTION  =  ■ 

THE  REDUCED  COSTS 

THE 

1-TH 

REDUCED  COST  = 

0.13878D-16 

THE 

2-TH 

REDUCED  COST  = 

-O.OOOOOD  00 

THE 

3-TH 

REDUCED  COST  = 

-O.OOOOOD  00 

THE 

4-TH 

REDUCED  COST  = 

0.44409D-15 

THE 

5-TH 

REDUCED  COST  = 

-O.OOOOOD  00 

THE 

6-TH 

REDUCED  COST  = 

-O.OOOOOD  00 

THE 

7-TH 

REDUCED  COST  = 

-0.90756D  00 

THE 

8-TH 

REDUCED  COST  = 

-O.OOOOOD  00 

THE 

9-TH 

REDUCED  COST  = 

-0,924370-01 

THE 

10-TH 

REDUCED  COST  = 

-O.OOOOOD  00 

THE 

11-TH 

REDUCED  COST  = 

-0.10000D  01 

THE 

12-TH 

REDUCED  COST  = 

-O.OOOOOD  00 

THE 

1 3-TH 

REDUCED  COST  = 

-0.10000D  01 

THE 

14-TH 

REDUCED  COST  = 

-O.OOOOOD  00 

THE 

15-TH 

REDUCED  COST  = 

-0.10000D  01 

THE 

16-TH 

REDUCED  COST  = 

-O.OOOOOD  00 

THE 

17-TH 

REDUCED  COST  = 

-0.92437D-01 

THE 

18-TH 

REDUCED  COST  = 

-O.OOOOOD  00 

THE 

19-TH 

REDUCED  COST  = 

-0.90756D  00 

TENTATIVELY  THE 

2-TH  BASIC  VARIABLE  IS  LEAVING  ' 

THE  YRJ'S 

f  YRf  J) 

=  0  IF 

THE  J-TH  VARIABLE  IS  IH  THE  BASIS 

IT  CAN  BE  ZERO  OTHERWISE  TOO.) 

YRf 

1)  = 

O.OOOOOD  00 

YRf 

2)  = 

O.OOOOOD  00 

YRf 

3)  = 

O.OOOOOD  00 

YRf 

4)  = 

-0.10000D  01 

YRf 

5)  = 

O.OOOOOD  00 

YRf 

6)  = 

O.OOOOOD  00 

YRf 

7)  = 

0.280110  00 

YRf 

8)  = 

O.OOOOOD  00 

YRf 

9)  = 

-0.2B011D  00 

YRf 

10)  * 

O.OOOOOD  00 

YRf 

11)  * 

O.OOOOOD  00 

YRf  12)  * 

O.OOOOOD  00 

-0.58047B  01 


-30- 

YR(  13)  =  0.00000D  00 

YRI  14)  *  O.OOOOOD  00 

YR(  IS)  =  O.OOOOOD  00 

YR(  14)  =  O.OOOOOD  00 

YR(  17)  =  -0.2801  ID  00 

YR<  18)  =  O.OOOOOD  00 

YR(  19)  =  0.2801  ID  00 

THE  RATIO  OF  THE  NET  PRICES  TO  THE  NEGATIVE  YRJ'S 
THE  RATIO  OF  THE  4-TH  NET  PRICE  TO  YR(  4)  = 

THE  RATIO  OF  THE  9-TH  NET  PRICE  TO  YR<  9)  = 

THE  RATIO  OF  THE  17-TH  NET  PRICE  TO  YR<  17)  = 

THE  2-TH  VARIABLE  IS  LEAVING  THE  BASIS. 

THE  4-TH  VARIABLE  IS  ENTERING  THE  BASIS. 

THE  BASIC  VARIABLES  ARE  NON 

10)  4>  5»  12»  14r  4)  14)  8)  18)  3> 

THE  VALUES  OF  THE  BASIC  VARIABLES  ARE  NON 


0.44409D-15 
-0.33000D  00 
-0.33000D  00 


XBI 

1)  = 

0.41977D  01 

XB< 

2)  = 

0.34497D-01 

XBI 

3)  = 

0.20988D  01 

XB( 

4)  = 

0.24942D-15 

XBI 

5)  = 

-0.70683D-15 

XBI 

4)  = 

0.32089D  01 

XBI 

7)  = 

0.44178D  01 

XBI 

8)  = 

0.49497D  00 

XBI 

9)  = 

0.99394D  00 

XBI 

10)  = 

0. 189210  00 

THE  CURRENT  VALUE  OF  THE  OBJECTIVE  FUNCTION  = 
THE  REDUCED  COSTS 


-0.58047D  01 


THE  1-TH  REDUCED  COST  =  0.27754D-14 

THE  2-TH  REDUCED  COST  =  -0.22204D-15 

THE  3-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  4-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  5-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  4-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  7-TH  REDUCED  COST  =  -0.90754D  00 

THE  8-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  9-TH  REDUCED  COST  =  -0.92437D-01 

THE  10-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  11-TH  REDUCED  COST  =  -0.10000B  01 

THE  12-TH  REDUCED  COST  *  -O.OOOOOD  00 

THE  13-TH  REDUCED  COST  =  -0.10000D  01 

THE  14-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  15-TH  REDUCED  COST  =  -0.10000D  01 

THE  14-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  17-TH  REDUCED  COST  =  -0.92437D-01 

THE  18-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  19-TH  REDUCED  COST  =  -0.90754D  00 

THE  CURRENT  BASIC  SOLUTION  IS  FEASIBLE  AND  HENCE  OPTIHAl. 

THE  NONZERO  VARIABLES  ARE  AS  FOLLOWS! 

X(  10)  =  0.41977D  01 

X(  4)  =  0.34497D-01 

X<  5)  =  0.20988D  01 

X(  12)  =  0.24942D-15 

X(  14)  =  -0.70683D-15 

XI  4)  =  0.32089D  01 

X(  14)  =  0.44178D  01 

X<  8)  =  0.49497D  00 

X<  18)  =  0.99394D  00 

XI  3)  =  0.18921D  00 

THE  OPTIMAL  VALUE  OF  THE  OBJECTIVE  FUNCTION  IS  -O.58047D  01 

SUPPLEMENTAL  INFORMATION  FROM  THE  MINI-MONTE  CARLO  STUDY 


THE  GENERATED  N(O)l)  EPSILONS 
-0.37319 
-0.57141 
0.04588 
0.10940 
0,29791 

THE  PROBLEM  P2  HAS  BEEN  SELECTED 
THE  INVERSE  OF  XTX 

0.810298  -0.148497 

-0.148497  0.044521 

VALUE  OF  BET AO  TO  COMPUTE  RHS 
-0.594250  0.134933 

THE  RHS)  YMXB)  FOR  PI  OR  P2 

0.047494 
-0.244174 
0.224595 
-0.039891 
-0.010224 


SAMPLE  NUMBER  = 


-31 


THE  INITIAL  VALUES  OF  THE  BASIC  VARIABLES 


XB< 

1)  = 

-0.56988D  00 

XB< 

2)  = 

-0.66336D-01 

XB( 

3)  = 

-0.28494D  00 

XB< 

4)  * 

0.802790-14 

XB( 

5)  = 

0.48323D-16 

XB( 

6)  = 

0.54175D  00 

XB( 

7)  = 

0.10835D  01 

XB( 

8)  = 

0.10662D  00 

XB( 

9)  * 

0.21323D  00 

XB( 

10)  * 

0.42693D  00 

THE  INITIAL  VALUE  OF  THE  OBJECTIVE  FUNCTION 
THE  REDUCED  COSTS 

THE  l-TH  REDUCED  COST 


-0.3634 3D  00 


THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 


2- TH  REDUCED  COST 

3- TH  REDUCED  COST 

4- TH  REDUCED  COST 

5- TH  REDUCED  COST 

6- TH  REDUCED  COST 

7- TH  REDUCED  COST 

8- TH  REDUCED  COST 

9- TH  REDUCED  COST 

10- TH  REDUCED  COST 
1  l-TH  REDUCED  COST 

12- TH  REDUCED  COST 

13- TH  REDUCED  COST 

14- TH  REDUCED  COST 

15- TH  REDUCED  COST 

16- TH  REDUCED  COST 

17- TH  REDUCED  COST 

18- TH  REDUCED  COST 

19- TH  REDUCED  COST 


0.27754P-14 
-0.22204D-15 
-O.OOOOOD  00 
-O.OOOOOD  00 
-O.OOOOOD  00 
-O.OOOOOD  00 
-0.90756D  00 
-O.OOOOOD  00 
-0.92437D-01 
-O.OOOOOD  00 
-0.10000D  01 
-O.OOOOOD  00 
-0.10000D  01 
-O.OOOOOD  00 
-0.10000D  01 
-O.OOOOOD  00 
— 0.92437D— 01 
-O.OOOOOD  00 
-0.90756D  00 


l-TH  BASIC  VARIABLE  IS  LEAVING  THE  BASIS 


THE 

TENTATIVELY  THE 
THE  YRJ'S 

(  YR(J)  =  0  IF  THE  J-TH  VARIABLE  IS  IN  THE  BASIS. 
IT  CAN  BE  ZERO  OTHERWISE  TOO. ) 


YR< 

YR< 

YR( 

YR( 

YR( 

YR( 

YR( 

YR( 

YR( 


1) 

2) 

3) 

4) 

5) 

6) 

7) 

8) 
9) 


YR(  10) 
YR<  11) 


YR(  12) 
YR(  13) 


YR(  14) 
YR(  15) 


YR( 

YR( 

YR( 

YR( 


16) 

17) 

18) 
19) 


O.OOOOOD  00 
0.22204D-15 
O.OOOOOD  00 
O.OOOOOD  00 

o.ooom  oo 

O.OOOOOD  00 
-0.30644D  01 
O.OOOOOD  00 
0.10644D  01 
O.OOOOOD  00 
O.OOOOOD  00 
O.OOOOOD  00 
O.OOOOOD  00 
O.OOOOOD  00 
-0.10000D  01 
O.OOOOOD  00 
0.30644D  01 
O.OOOOOD  00 
-0.10644D  01 


THE  RATIO  OF  THE  NET  PRICES  TO  THE  NEGATIVE  YRJ'S 


THE  RATIO  OF  THE  7-TH  NET  PRICE  TO  YR( 

THE  RATIO  OF  THE  15-TH  NET  PRICE  TO  YR( 

THE  RATIO  OF  THE  19-TH  NET  PRICE  TO  YR< 

l-TH  VARIABLE  IS  LEAVING  THE  BASIS. 

7-TH  VARIABLE  IS  ENTERING  THE  BASIS. 

THE  BASIC  VARIABLES  ARE  NOW 

7>  4i  5*  12»  14*  6,  16,  8, 

THE  VALUES  OF  THE  BASIC  VARIABLES  ARE  NOW 


7)  = 
15)  = 
19)  = 


THE 

THE 


-0.29616D  00 
-0.10000D  01 
-0.85263D  00 


18,  3, 


XB( 

1)  = 

0.18596D  00 

XB( 

2)  = 

-0.14245D-01 

XB( 

3)  = 

0.27756D-16 

XB( 

4)  = 

0.37193P  00 

XB( 

5)  = 

0.86648D-17 

XB( 

6)  = 

0.30005D  00 

XB( 

7)  = 

0.60009D  00 

XB( 

8)  = 

0.46190D-01 

XB< 

9)  = 

0.92380D-01 

XB( 

10)  = 

0.83651D-01 

THE  CURRENT  VALUE  OF  THE  OBJECTIVE  FUNCTION  * 
THE  REDUCED  COSTS 

THE  l-TH  REDUCED  COST  = 

“  2-TH  REDUCED  COST  = 

3- TH  REPUCED  COST  * 

4- TH  REDUCED  COST  = 

5- TH  REPUCED  COST  * 


-0.53220D  00 


THE 

THE 

THE 

THE 


-0.36082D-15 
-0.22204D-15 
-O.OOOOOD  00 
-O.OOOOOD  00 
-O.OOOOOD  00 


-32- 


THE  6-TH  REDUCED  COST  = 
THE  7-TH  REDUCED  COST  = 
THE  8-TH  REDUCED  COST  = 
THE  9-TH  REDUCED  COST  = 


THE  10-TH  REDUCED  COST  = 


THE  11-TH  REDUCED  COST  = 


THE  12-TH  REDUCED  COST  = 


-O.OOOOOD  00 
-O.OOOOOD  00 
-O.OOOOOD  00 
-0.40768D  00 
-0.296160  00 
-0.10000D  01 
-O.OOOOOD  00 


THE  13-TH  REDUCED  COST  =  -0.10000D  01 

THE  14-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  15-TH  REDUCED  COST  =  -0.70384D  00 

THE  16-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  17-TH  REDUCED  COST  =  -0.10000D  01 

THE  18-TH  REDUCED  COST  -  -O.OOOOOD  00 

THE  19-TH  REDUCED  COST  =  -0.59232D  00 

TENTATIVELY  THE  2-TH  BASIC  VARIABLE  IS  LEAVING  THE  BASIS 


THE  YRJ'S 

I  YR(J)  =  0  IF  THE  J-TH  VARIABLE  IS  IN  THE  BASIS. 


IT  CAN  BE  ZERO  OTHERWISE  TOO.) 

YR( 

1) 

= 

-0.13878D-16 

YR( 

2) 

= 

-0.10000D  01 

YR( 

3) 

= 

O.OOOOOD  00 

YR( 

4) 

= 

O.OOOOOD  00 

YR( 

5) 

= 

O.OOOOOD  00 

YR( 

6) 

- 

O.OOOOOD  00 

YR( 

7) 

S 

O.OOOOOD  00 

YR( 

8) 

= 

O.OOOOOD  00 

YRI 

9) 

r 

0.18282D  00 

YRI 

10) 

= 

-0.91408D-01 

YR( 

11) 

s 

O.OOOOOD  00 

YRI 

12) 

2 

O.OOOOOD  00 

YRI 

13) 

= 

O.OOOOOD  00 

YRI  14) 

r 

O.OOOOOD  OO 

YRI 

15) 

= 

0.91408D-01 

YRI 

16) 

s 

O.OOOOOD  00 

YRI 

17) 

= 

O.OOOOOD  00 

YRI 

18) 

= 

O.OOOOOD  00 

YRI 

19) 

s 

-0.18282D  00 

THE  RATIO  OF  THE  NET  PRICES  TO  THE  NEGATIVE  YRJ'S 

THE  RATIO  OF  THE  2-TH  NET  PRICE  TO  YK<  2)  =  -0.22204D-15 

THE  RATIO  OF  THE  10-TH  NET  PRICE  TO  YR(  10)  =  -0.32400D  01 

THE  RATIO  OF  THE  19-TH  NET  PRICE  TO  YR(  19)  =  -0.32400D  01 

THE  2-TH  VARIABLE  IS  LEAVING  THE  BASIS. 

THE  2-TH  VARIABLE  IS  ENTERING  THE  BASIS. 


THE  BASIC  VARIABLES  ARE  NOW 

7,  2,  5.  12.  14.  6.  16,  8,  18.  3. 

THE  VALUES  OF  THE  BASIC  VARIABLES  ARE  NOW 
XB(  1)  =  0.1B596D  00 

XB(  2)  =  0.14245D-01 

XB(  3)  =  0.38305D-16 

XB(  4)  =  0.37193D  00 

XB(  5)  =  0.18154D-16 

XB<  6)  =  0.30005D  00 

XB(  7)  =  0.60009D  00 

XBI  8)  =  0.46190D-01 

XB(  9)  =  0.92380D-01 

XB(  10)  =  0.83651D-01 

THE  CURRENT  VALUE  OF  THE  OBJECTIVE  FUNCTION  =  -0.53220D  00 

THE  REDUCED  COSTS 

THE  1-TH  REDUCED  COST  =  -0.31919D-15 


THE  2-TH  REDUCED  COST  *  -O.OOOOOD  00 

THE  3-TH  REDDED  COST  =  -O.OOOOOD  00 

THE  4-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  5-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  6-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  7-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  8-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  9-TH  REDUCED  COST  =  -0.40768D  00 

THE  10-TH  REDUCED  COST  =  -0.29616D  00 

THE  11-TH  REDUCED  COST  =  -0.10000D  01 

THE  12-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  13-TH  REDUCED  COST  =  -0.10000D  01 

THE  14-TH  REDUCED  COST  *  -O.OOOOOD  00 

THE  15-TH  REDUCED  COST  *  -0.70384D  00 

THE  16-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  17-TH  REDUCED  COST  =  -0.10000D  01 

THE  18-TH  REDUCED  COST  *  -O.OOOOOD  00 

THE  19-TH  REDUCED  COST  *  -0.59232D  00 

THE  CURRENT  BASIC  SOLUTION  IS  FEASIBLE  AND  HENCE  OPTIHAl, 
THE  NONZERO  VARIABLES  ARE  AS  FOLLOWS! 

XI  7)  =  0.18596D  00 

X<  2)  *  0.14245D-01 


X<  5)  =  0.38305D-16 

XI  12)  =  0.37193D  00 

X(  14)  =  0.18154B-16 

X(  6)  =  0.30005D  00 

X<  16)  =  0.60009D  00 

XI  8)  =  0.46190D-01 

X(  18)  =  0.92380D-01 

X(  3)  =  0.83651D-01 

THE  OPTIMAL  VALUE  OF  THE  OBJECTIVE  FUNCTION  IS 


-0.53220D  00 


SUPPLEMENTAL  INFORMATION  FROM  THE  MINI-MONTE  CARLO  STUDY 


THE  GENERATED  NIO.l)  EPSILONS  SAMPLE  NUMBER  =  2 

0,87619 
0,59391 
-0.35016 
0.51585 
1.51357 

THE  PROBLEM  PI  HAS  BEEN  SELECTED 
THE  INVERSE  OF  XTX 

0.810298  -0.168497 

-0.168497  0.046521 

VALUE  OF  BET AO  TO  COMPUTE  RHS 
0.216974  0.113998 

THE  RHSr  YMXB.  FOR  PI  DR  P2 
-0.531540 
-0.154642 
0.911404 
0.320132 
-0.545353 

THE  INITIAL  VALUES  OF  THE  BASIC  VARIABLES 
XB<  1)  =  0.14477D  01 

XB<  2)  =  0.25252D-02 

XB<  3)  =  -0.19987D-15 

XB(  4)  =  0.28955D  01 

XB(  5)  =  — 0. 13871D-15 

XB(  6)  =  -0.37899D  00 

XB(  7)  =  -0.75799D  00 

XB(  8)  =  -0.86256D  00 

XB(  9)  =  -0.17251D  01 

XB(  10)  =  -0.52871D  00 

THE  INITIAL  VALUE  OF  THE  OBJECTIVE  FUNCTION  =  -0.20619B  00 

THE  REDUCED  COSTS 

THE  1-TH  REDUCED  COST  =  -0.31919D-15 

THE  2-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  3-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  4-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  5-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  6-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  7-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  8-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  9-TH  REDUCED  COST  =  -0.40768D  00 

THE  10-TH  REDUCED  COST  =  -0.29616D  00 

THE  1 1-TH  REDUCED  COST  =  -0.10000D  01 

THE  12-TH  REDUCED  COST  =  -O.OOOOOD  00 


THE  13-TH  REDUCED  COST  =  -0.10000D  01 

THE  14-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  15— TH  REDUCED  COST  =  -0.70384D  00 

THE  16-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  17-TH  REDUCED  COST  =  -0.10000D  01 

THE  18-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  19-TH  REDUCED  COST  =  -0.59232D  00 

TENTATIVELY  THE  9-TH  BASIC  VARIABLE  IS  LEAVING  THE  BASIS 
THE  YRJ'S 

(  YR(J)  =  0  IF  THE  J-TH  VARIABLE  IS  IN  THE  BASIS. 

IT  CAN  BE  ZERO  OTHERWISE  TOO.) 


YRI  10)  = 
YR(  11)  = 
YRI  12)  = 
YRI  13)  * 
YRI  14)  * 
YRI  15)  = 


0.44409D-15 
O.OOOOOD  00 
O.OOOOOD  00 
0.22204D-15 
O.OOOOOD  00 
O.OOOOOD  00 
O.OOOOOD  00 
O.OOOOOD  00 
0.15759D  01 
0.21207D  00 
O.OOOOOD  00 
O.OOOOOD  00 
-O.IOOOOD  01 
O.OOOOOD  00 
-0.21207D  00 


(etc. 


THE  VALUES  OF  THE  BASIC  VARIABLES  ARt  NON 
XB(  1)  =  -0.41633D-16 

XB(  2)  =  0.56469D-01 

XB(  3)  =  0.73958P  00 

XB<  4)  «  0.34616D  01 

X8(  5)  =  0.18851D  00 

XB(  6)  =  0 .173081*  01 

XB<  7)  =  0.249200-14 

XB(  8)  =  0.942540-01 

XB(  9)  =  0.147920  01 

XB(  10)  =  0.383920  00 

THE  CURRENT  VALUE  OF  THE  OBJECTIVE  FUNCT.ON  = 
THE  REDUCED  COSTS 


-0.254440  01 


THE  1-TH  REDUCED  COST  =  -0.112410-14 

THE  2-TH  REDUCED  COST  -  -0.000000  00 

THE  3-TH  REDUCED  COST  =  -0.000000  00 

THE  4-TH  REDUCED  COST  =  -0,444090-15 

THE  5-TH  REDUCED  COST  =  -0.000000  00 

THE  4-TH  REDUCED  COST  =  -0.000000  00 

THE  7-TH  REDUCED  COST  =  -0.000000  00 

THE  8-TH  REDUCED  COST  =  -0.000000  00 

THE  9-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  10-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  11-TH  REDUCED  COST  =  -0.10000D  01 

THE  12-TH  REDUCED  COST  =  -0.77871D  00 

THE  13-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  14-TH  REDUCED  COST  =  -0.72129D  00 

THE  15-TH  REDUCED  COST  =  -0.10000P  01 

THE  14-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  17-TH  REDUCED  COST  =  -0.22129D  00 

THE  18-TH  REDUCED  COST  =  -0.10000D  01 

THE  19-TH  REDUCED  COST  =  -0.27871B  00 

THE  CURRENT  BASIC  SOLUTION  IS  FEASIBLE  AND  HENCE  OPTIMAL. 

THE  NONZERO  VARIABLES  ARE  AS  FOLLOWS: 

X<  7)  =  -0.41433B-14 

X(  2)  =  0.54449D-01 

X(  5)  =  0.73958D  00 

X<  14)  =  0.34414D  01 

X(  13)  =  0.18851D  00 

X(  4)  =  0.17308D  01 

X<  9)  =  0.24920D-16 

X<  8)  =  0. 94254D-01 

X(  10)  =  0.14792D  01 

X(  3)  =  0.38392D  00 

THE  OPTIMAL  VALUE  OF  THE  OBJECTIVE  FUNCTION  IS  -0.25444 


-0.25446B  01 


SUPPLEMENTAL  INFORMATION  FROM  THE  MINI-MONTE  CARLO  STUDY 


THE  GENERATED  N(0,1)  EPSILONS  SAMPLE  NUMBER  = 
0.41901 
-0.48994 
-1.15302 
-0.01590 
-1.46441 

THE  PROBLEM  P2  HAS  BEEN  SELECTED 
THE  INVERSE  OF  XTX 

0.ei0298  -0.168497 

-0.168497  0.044521 

VALUE  OF  BETAO  TO  COMPUTE  RHS 
0.178028  -0.209520 

THE  RHSf  TMXBt  FOR  PI  OR  P2 
0.475445 
-0.259401 
-0.698296 
0.943740 
-0.461707 

THE  INITIAL  VALUES  OF  THE  BASIC  VARIABLES 
XB(  l)  *  -0.17131D-15 

XB(  2)  »  -C.66271D-01 

XB<  3)  *  0.12999D  01 

XB<  4)  «  -0.10196D  01 

X8<  5)  *  0.29647D  01 

xf!  $1  *  ^.'ll43f-?6 

XB(  8)  =  0.14823D  01 

XB(  9)  *  0.25997D  01 

XB(  10)  -  -0.B9844D  00 

THE  INI  UAL  VALUE  OF  THE  OBJECTIVE  FUNCTION  *  - 

THE  REDUCED  COSTS 

THE  1-TH  REDUCED  COST  *  -0.11241D-14 
THE  2-TH  REDUCED  COST  =  -O.OOOOOD  00 


-0*227240  01 
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THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 


3- TH  REDUCED  COST 

4- TH  REDUCED  COST 

5- TH  REDUCED  COST 

6- TH  REDUCED  COST 

7- TH  REDUCFD  COST 

8- TH  REDUCED  COST 

9- TH  REDUCED  COST 

10- TH  REDUCED  COST 

11- TH  REDUCED  COST 

12- TH  REDUCED  COST 

13- TH  REDUCED  COST 

14- TH  REDUCED  COST 

15- TH  REDUCED  COST 
14-TH  REDUCED  COST 

17- TH  REDUCED  COST 

18- TH  REDUCED  COST 

19- TH  REDUCED  COST 


THE 

TENTATIVELY  THE 
THE  YRJ'S 

(  YR(J)  =  0  IF  THE  J-TH  VARIABLE 


-O.OOOOOD  00 
-0.44409D-15 
-O.OOOOOD  00 
-O.OOOOOD  00 
-O.OOOOOD  00 
-O.OOOOOD  00 
-O.OOOOOD  00 
-O.OOOOOD  00 
-0.10000D  01 
-0.77871D  00 
-O.OOOOOD  00 
-0.72129D  00 
-0.10000D  01 
-O.OOOOOD  00 
-0.22129D  00 
-0.10000D  01 
-0.27871D  00 


4-TH  BASIC  VARIABLE  IS  LEAVING  THE  BASIS 


IS  IN  THE  BASIS. 


CAN  BE  ZERO 

OTHERWISE  TOO.) 

YR( 

1) 

= 

0.51348D-15 

YR( 

2) 

S 

O.OOOOOD  00 

YR< 

3) 

= 

O.OOOOOD  00 

YR( 

4) 

s 

-0.11102D-14 

YR( 

5) 

s 

O.OOOOOD  00 

YR( 

4) 

= 

O.OOOOOD  00 

YR( 

7) 

s 

O.OOOOOD  00 

YR< 

8) 

s 

O.OOOOOD  00 

YR( 

9) 

= 

O.OOOOOD  00 

YR< 

10) 

r 

O.OOOOOD  00 

YR( 

11) 

= 

-0.10000D  01 

YR( 

12) 

r 

0.12997D  01 

YR< 

13) 

= 

O.OOOOOD  00 

YR< 

14) 

= 

-0.29972D  00 

YR< 

15) 

= 

O.OOOOOD  00 

YR( 

14) 

O.OOOOOD  00 

YR( 

17) 

= 

-0.12997D  01 

YR( 

18) 

= 

O.OOOOOD  00 

YR( 

19) 

= 

0.29972D  00 

:  RATIO  OF 

THE  NET  PRICES  TO 

THE  RATIO  OF  THE  11-TH  NET  F'RICE  TO  YR( 

THE  RATIO  OF  THE  14-TH  NET  FRICE  TO  YR< 

THE  RATIO  OF  THE  17-TH  NET  FRICE  TO  YR( 

4-TH  VARIABLE  IS  LEAVING  THE  BASIS. 
17-TH  VARIABLE  IS  ENTERING  THE  BASIS. 

THE  BASIC  VARIABLES  ARE  NOW 

7,  2.  5i  17,  13,  4,  9,  8, 

THE  VALUES  OF  THE  BASIC  VARIABLES  ARE  NOW 


THE 

THE 


11)  S 
14)  * 
17)  = 


-0.10000D  01 
-0.24045D  01 
-0.17026D  00 


10,  3, 


XB( 

1)  = 

0.39224D  00 

XB( 

2)  = 

0.43400D-01 

XB( 

3)  = 

0.49884D  00 

X8( 

4)  = 

0.78448D  00 

XB( 

5)  = 

0.27098D  01 

XB( 

4)  = 

-0.11102D-15 

XB( 

7)  = 

0.72940D-14 

XB( 

8)  = 

0.13549D  01 

XB< 

9)  = 

0.13977B  01 

XB( 

10)  = 

-0.17438D  00 

CURRENT  VALUE  OF  THE  OBJECTIVE  FUNCTION  = 

THE  REDUCED  COSTS 


-0.24440D  01 


THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 


W 


THE 

THE 

THE 


1- TH  REDUCED  COST  = 

2- TH  REDUCED  COST  = 

3- TH  REDUCED  COST  = 

4- TH  REDUCED  COST  = 

5- TH  REDUCED  COST  = 
4-TH  REDUCED  COST  = 

7- TH  REDUCED  COST  = 

8- TH  REDUCED  COST  = 

9- TH  REDUCED  COST  = 

10- TH  REDUCED  COST  = 

11- TH  REDUCED  COST  = 

12- TH  REDUCED  COST  = 

13- TH  REDUCED  COST  = 

14- TH  REDUCED  COST  * 

15- TH  REDUCED  COST  = 
14-TH  REDUCED  COST  = 

17- TH  REDUCED  COST  = 

18- TH  REDUCED  COST  = 

19- TH  REDUCED  COST  = 


-0.99920D-15 
-O.OOOOOD  00 
-O.OOOOOD  00 
-O.OOOOOD  00 
-O.OOOOOD  00 
-O.OOOOOD  00 
-O.OOOOOD  00 
-O.OOOOOD  00 
-O.OOOOOD  00 
-O.OOOOOD  00 
-0.82974D  00 
-0.10000D  01 
-O.OOOOOD  00 
-0.47024D  00 
-0.10000D  01 
-0.17024D  00 
-O.OOOOOD  00 
-0.10000D  01 
-0.32974D  00 


TENTATIVELY  THE  10-TH  BASIC  VARIABLE  IS  LEAVING  THE  BASIS 


'L  . 
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THE  YRJ'S 
(  YR(J)  =  0  IF  THE 


IT  CAN  BE  ZERO  OTHERWISE  TOO.) 

YR( 

1) 

s 

-0.10000D  01 

YR( 

2) 

= 

O.OOOOOD  00 

YR( 

3) 

= 

O.OOOOOD  00 

YR( 

4) 

= 

0.66613D-15 

YR( 

5) 

s 

O.OOOOOD  00 

YR( 

6) 

2 

O.OOOOOH  00 

YR( 

7) 

= 

O.OOOOOD  00 

YR( 

8) 

= 

O.OOOOOD  00 

YR( 

9) 

= 

O.OOOOOD  00 

YR( 

10) 

S 

O.OOOOOD  00 

YR< 

11) 

2 

0 . 7101 3D  00 

YR< 

12) 

= 

0.18041D-15 

YR( 

13) 

= 

O.OOOOOD  00 

YR( 

14) 

Z 

-0,21 01 3D  00 

YR( 

15) 

= 

O.OOOOOD  00 

YR( 

16) 

z 

-0.71013D  00 

YR< 

17) 

= 

O.OOOOOD  00 

YR( 

18) 

z 

O.OOOOOD  00 

YR( 

19) 

z 

0.21013D  00 

THE  RATIO  OF 

THE  NET  PRICES  TO 

THE 

THE 


THE  RATIO  OF  THE  1-TH  NET  PRICE  TO  YR<  1)  =  -0.99920D-15 
THE  RATIO  OF  THE  14-TH  NET  PRICE  TO  YR(  14)  =  -0.31S97D  01 
THE  RATIO  OF  THE  16-TH  NET  PRICE  TO  YR(  16)  =  -0.23976D  00 
10— TH  VARIABLE  IS  LEAVING  THE  BASIS. 

1-TH  VARIABLE  IS  ENTERING  THE  BASIS. 


THE  BASIC  VARIABLES  ARE  NON 


7) 

2) 

5)  17)  13) 

6>  9>  8>  10) 

VALUES  OF  THE  BASIC  VARIABLES  ARE  HOW 

XB( 

1)  = 

0.39224D  00 

XB( 

2)  = 

0.43600D-01 

XB( 

3)  = 

0.69886D  00 

XB( 

4)  = 

0.7B448D  00 

XB( 

5)  = 

0.27098D  01 

XB( 

6)  = 

-0.12796D-15 

XB( 

7)  = 

0.85041D-16 

XB( 

8)  = 

0.13549D  01 

XB( 

9)  = 

0.13977D  01 

XB(  10)  = 

0.17438D  00 

CURRENT  VALUE  OF  THE  OBJECTIVE  FUNCTION  = 

REDUCED  COSTS 

THE 

1-TH 

REDUCED  COST  = 

-O.OOOOOD  00 

THE 

2-TH 

REDUCED  COST  = 

-O.OOOOOD  00 

THE 

3-TH 

REDUCED  COST  = 

-0.12490D-15 

THE 

4-TH 

REDUCED  COST  = 

-O.OOOOOD  00 

THE 

5-TH 

REDUCED  COST  = 

-O.OOOOOD  00 

THE 

6-TH 

REDUCED  COST  = 

-O.OOOOOD  00 

THE 

7-TH 

REDUCED  COST  = 

-O.OOOOOD  00 

THE 

8-TH 

REDUCED  COST  = 

-O.OOOOOD  00 

THE 

9-TH 

REDUCED  COST  = 

-O.OOOOOD  00 

THE 

10-TH 

REDUCED  COST  = 

-O.OOOOOD  00 

THE 

11-TH 

REDUCED  COST  = 

-0.B2974D  00 

TIC 

12-TH 

REDUCED  COST  = 

-O.IOOOOD  01 

THE 

13-TH 

REDUCED  COST  = 

-O.OOOOOD  00 

THE 

14-TH 

REDUCED  COST  = 

-0.67026D  00 

THE 

15-TH 

REDUCED  COST  = 

-O.IOOOOD  01 

THE 

16-TH 

REDUCED  COST  = 

-0.17026D  00 

THE 

17-TH 

REDUCED  COST  = 

-O.OOOOOD  00 

THE 

18-TH 

REDUCED  COST  = 

-O.IOOOOD  01 

THE 

19-TH 

REDUCED  COST  = 

-0.32974D  00 

-0.24460B  01 


THE  CURRENT  BASIC  SOLUTION  IS  FEASIBLE  AND  HENCE  OPTIMAL. 
THE  NONZERO  VARIABLES  ARE  AS  FOLLOWS! 


X< 

7)  = 

0.39224D  00 

X( 

2)  = 

0.43600D-01 

X< 

5)  = 

0.69886D  00 

x< 

17)  = 

0.78448D  00 

X<  13)  = 

0.27098D  01 

X< 

6)  = 

-0.12796D-15 

X< 

9)  = 

0.85041B-16 

XI 

8)  = 

0.13549D  01 

X(  10)  = 

0.13977D  01 

X< 

1)  = 

0.17438D  00 

THE  OPTIKAL  VALUE  OF  THE  OBJECTIVE  FUNCTION  IS  -0.24460D  01 


THE  AUXILIARY  LEAST  SOUARES  ESTIMATE)  BETAOi  OF  THE  REGRESSION  PARAMETER  VECTOR)  BETA 
LEAST  SQUARES  ESTIMATE  OF  BETA!  1)  =  3.227002 

LEAST  SQUARES  ESTIMATE  OF  BETA!  2)  *  -1.159747 

MRS.  A'S  ANSWER  !  THE  ESTIMATE  OF  THE  REGRESSION  PARAMETER  VECTOR  WHICH  MINIMIZES  THE  SUM  OF  THE  ABSOLUTE  RESIBUALS! 


LI  ESTIMATE  OF  BETAt  1)  =  3.416213 

LI  ESTIMATE  OF  BETA!  2)  =  -1.12324? 

THE  RESIDUALS.  R(I).I=l»NOBS 
2.098826 
3.208877 
0.000000 
0.49696? 

0.000000 

THE  SUH  OF  THE  ABSOLUTE  RESIDUALS  =  5.804672 

THE  MAXIMUM  ABSOLUTE  RESIDUAL  =  3.208877 

AUXILIARY  RESULTS  OF  THE  MINI-MONTE  CARLO  STUDY 

VALUES  OF  DELTA  BETA  STAR 


SAMPLE  NUMBER 

= 

l 

-0.51060 

0.12269 

SAMPLE  NUMBER 

= 

2 

0.63765 

-0.02243 

SAMPLE  NUMBER 

S 

3 

0.94741 

-0.45843 

SAMPLE  NUMBER 

S 

4 

0.05573 

-0.14477 

SAMPLE  NUHBER 

X 

5 

-0.08930 

-0.16773 

SAMPLE  NUMBER 

= 

6 

-1.10362 

0.22395 

SAMPLE  NUHBER 

s 

7 

-0.86781 

-0.04936 

SAMPLE  NUMBER 

s 

8 

-1.24035 

0.54421 

SAMPLE  NUHBER 

= 

9 

0.44742 

0.24662 

SAMPLE  NUHBER 

= 

10 

-1.03640 

0.36600 

SAMPLE  NUMBER 

= 

11 

0.83504 

0.01230 

SAMPLE  NUMBER 

= 

12 

0.32935 

-0.04047 

SAMPLE  NUMBER 

= 

13 

-1.27719 

0.27656 

SAMPLE  NUMBER 

= 

14 

-0.58869 

0.33108 

SAMPLE  NUMBER 

= 

15 

-0.04885 

0.14684 

SAMPLE  NUMPER 

= 

16 

0.98576 

-0.03155 

SAMPLE  NUMBER 

s 

17 

-1.01343 

0.04536 

SAMPLE  NUMBER 

= 

18 

1.70310 

-0.28731 

SAMPLE  NUMBER 

= 

19 

-0.37732 

-0.04003 

SAMPLE  NUMBER 

= 

20 

0.00365 

-0.25312 

ESTIMATED  COVARIANCE  OF  DELTA  BETA  STAR 
0.774642  -0.148480 

-0.148480  0.063634 

SUH  OF  THE  OPTIMAL  OBJECTIVE  FUNCTIONS  OVER  ALL  SAMPLES  = 


48,526882 


AUXILIARY  RESULT*.  SIGMA  HAT  5  =  2.425026 

MAIN  RESULTS  OF  THE  HINI-MONTE  CARLO  STUDY 
ESTIMATED  VALUE  OF  SIGMA  (SIGMA  HAT  4)  =  2.392353 

ESTIMATED  COVARIANCE  OF  THE  REGRESSION  PARAMETER  VECTOR  (BETA)  USING  THIS  ESTIMATE  OF  SIGH,' 
4.433552  -0.849803 
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